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LIST  OF  SYMBOLS 


M  =   Random  variable  of  the  number  of  survivors 

for  the  X  force  in  a  stochastic  model. 

N  =   Random  variable  of  the  number  of  survivors 

for  the  Y  force  in  a  stochastic  model. 

=   Number  of  survivors  of  X  force  at  time  t 
for  a  deterministic  model. 

=   Number  of  survivors  of  Y  force  at  time  t 
for  a  deterministic  model. 

=  Realization  of  M  at  time  t. 

=  Realization  of  N  at  time  t. 

=  Initial  size  of  X  force. 

=  Initial  size  of  Y  force. 

=   Ratio  of  the  number  of  casualties  to  the 
initial  force  size  of  the  X  force  at  which 
the  X  force  disengages. 

FY  —   Ratio  of  the  number  of  casualties  to  the 

initial  force  size  of  the  Y  force  at  which 
the  Y  force  disengages. 

m,  =   Largest  integer  part  of  (1-F  ) *m  .   The  number 

p  of  survivors  at  which  the  X  force  disengages. 

n,  =   Largest  integer  part  of  (1-F„)  *n  .   The  number 

"  of  survivors  at  which  the  Y  force  disengages. 

F(t,x,y)     =   Rate  of  attrition  of  the  X  force  for  a 

deterministic  model. 

G(t,x,y)     =   Rate  of  attrition  of  the  Y  force  for  a 

deterministic  model. 

a  =   Rate  at  which  a  Y  combatant  attrits  the  X  force 

b  =   Rate  at  which  a  X  combatant  attrits  the  Y  force 


x(t) 

y(t) 

m(t) 

n(t) 

mo 

n 
o 

F 
X 
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F(t,m,n)dt   =   Probability  that  in  the  small  time  interval 

dt,  an  X  combatant  is  destroyed,  the  number 
of  fighting  units  existing  at  the  time  t 
preceeding  that  interval  being  m  and  n; 
for  a  stochastic  model. 

G(t,m,n)dt   =   Probability  that  in  the  small  time  interval 

dt,  a  Y  combatant  is  destroyed,  the  number 
of  fighting  units  existing  at  the  time  t 
preceeding  that  interval  being  m  and  n; 
for  a  stochastic  model. 

A(m,n)       =   Special  case  of  F(t,m,n)  with  constant 

coefficients . 

B(m,n)       =   Special  case  of  G(t,m,n)  with  constant 

coefficients . 

m(t)  =  Expected  value  of  M  at  time  t. 

n(t)  =  Expected  value  of  N  at  time  t. 

V(M,t)  =  Variance  of  M  at  time  t. 

V(N,t)  =  Variance  of  N  at  time  t. 

P(t,m,n)     =   Probability  that  at  time  t,  X  has  m  survivors 

and  Y  has  n  survivors,  given  that  at  time 

t  =  0 ,  m  =  m   and  n  =  n  . 
o  o 

P  =   Probability  X  wins  in  a  stochastic  model. 

A 

P  =  Probability  Y  wins  in  a  stochastic  model. 

P(m,n;r,s)   =   Probability  that  the  X  force  has  m  survivors 

and  the  Y  force  has  n  survivors  given  that 
they  started  with  r  combatants  and  s  combatants, 
respectively . 

Ax(t)       =  m(t)  -  x(t) 

Ay(t)        =  n(t)  -  y(t) 

A%y(t)       =  bias  for  the  X  forces  of  the  two 

models  at  time  t. 

A%y(t)       =   Percent  of  bias  for  the  Y  forces  of  the  two 

models  at  time  t. 
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I.   INTRODUCTION 

Although  combat  is  a  complex  random  process,  analysts 
frequently  model  it  with  deterministic  Lanchester-type 
equations  for  reasons  of  mathematical  tractability  and  com- 
putational convenience.   Such  deterministic  models  are  being 
used  extensively  today  in  defense  planning  studies.   However, 
one  can  include  the  various  random  aspects  of  combat  (e.g., 
uncertainty  in  enemy  force  level,  random  occurrence  of 
casualties,  etc.)  in  models  and  develop  so-called  stochastic 
models  of  the  combat  process.   Thus,  a  basic  way  of  classi- 
fying a  combat  model  is  whether  the  model  is  a  deterministic 
model  or  a  stochastic  one. 

Moreover,  there  are  fundamental  differences  between 
deterministic  and  stochastic  models.   Deterministic  dynamic 
combat  models  predict  the  future  with  certainty;  for  given 
initial  conditions  there  is  no  question  about  what  the  state 
of  the  conflict  will  be  at  any  future  time.   On  the  other 
hand,  stochastic  dynamic  models  only  tell  one  the  "chances" 
of  what  will  happen.   With  a  stochastic  model  one  does  not 
know  the  future  with  certainty.   That  is,  it  cannot  be 
guaranteed  that  a  specific  state  will  be  reached  in  the 
future  —  there  is  only  a  "probability"  it  will  be  reached. 

Each  type  of  model  has  inherent  advantages  and  disadvantages 
which  are  partially  summarized  in  Table  I. 
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Deterministic  Models 


Advantages 


Disadvantages 


1.  Requires  one  solution  of  the 
model  for  a  useful  outcome 

2.  Less  expensive  to  run 

3.  Interactions  are  more 
easily  analyzed 


1.   A  further  abstraction 
from  reality 


Stochastic  Models 


Advantages 


Disadvantages 


1.   Closer  to  reality 


1.  Requires  many  replica- 
tions to  provide 
statistically 
significant  results 

2.  More  expensive  to  run 

3.  Interactions  are  much 
more  difficult  to 
analyze 


TABLE  I.   Advantages  and  Disadvantages 

of  Deterministic  and  Stochastic 
Models 
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The  solutions  to  stochastic  models  quickly  become  quite  com- 
plex as  the  number  of  combatants  increases,  and  are  not  very 
useful;  so  one  must  use  other  methods  to  obtain  insights  into 
the  dynamics  of  combat.   The  two  most  commonly  used  methods, 
Monte  Carlo  simulation  and  finite  difference  approximation, 
both  possess  certain  disadvantages.   Monte  Carlo  combat  simu- 
lations require  that  the  battle  be  replicated  many  times  in 
order  to  obtain  a  good  statistical  estimate  of  the  probable 
course  of  battle,  while  finite  different  approximation  methods 
usually  require  digital  computer  implementation  (including 
development  of  a  computer  program) . 

On  the  other  hand,  solutions  to  deterministic  models  can 
sometimes  be  useful  for  developing  insights  into  the  dynamics 
of  combat.   Even  when  they  are  not,  the  model  must  be  solved 
only  once  by  finite  difference  approximation  as  opposed  to 
the  many  replications  required  for  Monte  Carlo  simulation. 
Also,  finite  difference  approximations  of  deterministic  models 
are  generally  less  costly  in  terms  of  time,  and  manhours 
required  to  implement  the  computer  program,  than  the  finite 
difference  approximations  of  stochastic  models.   The  economic 
costs  of  using  the  models  then,  is  the  basic  reason  most  models 
of  combat  used  in  the  past  have  been  deterministic.   But  can 
the  complex  random  process  of  combat  be  adequately  represented 
by  a  deterministic  model?   Does  one  destroy  any  of  the  essential 
features  of  the  random  combat  process  by  considering  a  deter- 
ministic model  as  representing  the  "average"  course  of  combat? 
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The  objective  of  this  thesis  is  to  consider  this  question  by 
investigating  the  differences  in  results  obtained  from  a 
deterministic  and  a  stochastic  model  of  combat. 

The  importance  of  this  objective  stems  from  the  fact  that 
both  types  of  models  are  used  extensively  today  in  the  decision 
making  process  in  the  U.S.  Army.   For  example,  the  Dynamic 
Tactical  Simulation,  DYNTACS-X,  and  the  Battalion  Differential 
Model,  BLDM,  widely  known  as  the  Bonder/IUA  model,  are  two 
models  that  simulate  battalion  level  mid-intensity  armored 
combat.   DYNTACS  is  an  event  sequenced  stochastic  model  that 
is  being  or  has  been  used  by  the  Rock  Island  Arsenal  to  evalu- 
ate proposed  mobility  improvements  to  the  M60  tank,  the  family 
of  scatterable  mines  concept,  cannon  launched  guided  projec- 
tiles, remotely  piloted  vehicles,  and  the  XM1  tank.   Bonder/ 
IUA  is  a  deterministic  Lanchester-type  model  that  is  being  or 
has  been  used  by  the  Rock  Island  Arsenal  in  sensitivity  analyses 
and  evaluation  of  the  MBT70  study,  several  anti-armor  automatic 
cannon  concepts,  and  the  Low  Dispersion  Automatic  Cannon  study. 
A  study  done  by  the  Rock  Island  Arsenal  [Ref.  5]  shows  that 
there  are  differences  in  the  results  of  the  two  models  with 
"equivalent"  inputs.   With  the  present  state  of  the  art,  it 
is  not  possible  to  predict  with  any  accuracy  when  significant 
differences  between  the  results  of  the  models  will  occur.   It 
is  felt  that  this  thesis  may  provide  trends  which  will  allow 
a  user  to  gain  insights  into  these  differences  and  more 
accurately  determine  when  a  deterministic  model  represents 
combat  "satisfactorily". 
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Models  have  had  wide  application  in  the  military  decision 
making  process,  and  the  realm  of  possible  future  applications 
is  even  larger.   For  example,  the  Engineer  Strategic  Studies 
Group  [Pef.  10]  has,  as  its  name  implies,  evaluated  several 
aspects  of  military  engineering,  but  with  little  emphasis  on 
the  use  of  combat  models.   Much  of  this  has  been  caused  by  a 
dearth  of  models  that  could  be  applied  to  engineer  studies. 
However,  this  is  being  altered;  for  example,  the  U.S.  Army 
Concepts  Analysis  Agency  [Ref.  30]  has  accepted  as  a  future 
tool  a  model  developed  by  ESSG,  the  Fortification  and  Obstacle 
Effect  Simulation,  1975.   This  is  a  high  resolution  model 
used  in  analyzing  obstacle  effectiveness.   The  model  is  part 
of  a  research  effort  by  ESSG  to  design  a  method  for  assessing 
the  contribution  that  obstacles  make  to  combat  operations. 
This  model  was  developed  using  ideas  from  other  models,  but 
it  was  based  on  the  identification  of  specific  objectives. 
From  these  objectives  measures  of  effectiveness  were  developed, 
and  from  these  measures  the  model  was  developed.   It  is  possi- 
ble, with  the  identification  of  proper  objectives,  to  develop 
models  for  other  phases  of  military  engineering  and,  more 
generally,  for  the  evaluation  of  other  phases  of  combat.   It 
is  therefore  felt  that  combat  models  do,  and  will  continue  to 
have  an  ever  increasing  application  to  the  evaluation  of 
military  systems. 

As  was  noted,  however,  it  is  not  always  clear  which  type 
of  model  should  be  applied  —  a  stochastic  or  a  deterministic 
model.   From  an  economic  standpoint,  the  advantages  of  the 
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deterministic  model  weight  heavily  in  its  favor.   If  the 
deterministic  results  do  not  significantly  differ  from  those 
of  the  stochastic  model,  it  seems  the  former  should  be  used. 
The  question  is,  when  are  the  models  "similar"? 
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II.   LITERATURE  SURVEY 

This  chapter  briefly  reviews  Lanchester's  well-known 
equations  of  combat  and  some  simple  extensions,  with  the 
discussion  focused  on  stochastic  formulations  of  the  combat 
attrition  process.   The  review  consists  of  three  parts; 
(A)  Deterministic  Models,  (B)  Stochastic  Models,  and 
(C)  Comparisons  of  the  two  types. 

A.   DETERMINISTIC  MODELS 

3 

In  1914  Lanchester  [Ref.  19]  hypothesized  that  under 
conditions  of  "modern  warfare"  combat  between  two  homogeneous 
forces  could  be  described  by  the  equations 


&  -  -  <*  «« 


$  -  -  bx  (2) 


with  initial  conditions 


x(t  =  0)  =  xQ  (3) 


y(t  =  0)  =  yQ  (4) 


where  a  and  b  are  commonly  referred  to  as  the  Lanchester 
attrition  rate  coefficients  and  x(t)  and  y(t)  are  force 
levels.   One  set  of  assumptions  that  has  been  hypothesized 
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6 

(proposed  by  H.  Weiss  [Ref.  35]  and  cited  by  Dolansky 
[Ref.  9])  to  yield  equations  (1)  and  (2)  is 

1.  Two  forces  are  engaged  in  a  battle.  The  units 
are  homogeneous ,  but  the  rate  of  attrition  may 
be  different  for  each  force. 

2.  Each  firer  on  either  side  is  within  the  range 
of  all  targets  on  the  other  side. 

3.  Attrition  rate  coefficients  are  constant. 

4(a).   Firers  have  perfect  knowledge  of  the  target 
locations  and  disposition  so  that  they  fire  only 
at  live  targets,  and  fire  can  be  shifted 
instantaneously. 

5(a).   Fire  is  uniformly  distributed  over  remaining 
targets . 

The  state  solution  (time  "independent"  solution  relating 

the  force  levels  )  is  given  by 


b(xQ2  -  x(t)2)  =  a(yQ2  -  y(t)2)  (5) 


hence  the  name  "square-law"  attrition  process. 

The  time  solution  for  equations  (1)  and  (2)  is  given 
by  (see  Morse  and  Kimball  [Ref.  22]) 


(t)  =  x  cosh  {-y   ab  t)  -   -\   a/b  y   sinh  (y  ab  t)  (6) 


y(t)  =  y   cosh  (\  ab  t)  -  \  b/a  x   sinh  (~\j   ab  t)  (7) 

The  square-law  was  developed  as  a  by-product  of  the 
original  work  by  Lanchester  to  quantitatively  justify  the 
principle  of  concentration  under  certain  combat  conditions 
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Lanchester  was  contrasting  these  conditions  with  the  condi- 
tions of  "ancient  warfare"  between  two  homogeneous  forces. 
He  postulated  that  the  latter  could  be  described  by  the 
equations 


at  =  -  axy  (8) 


§£  =  -  bxy  (9) 


with  initial  conditions 


x(t  =  0)  =  x  (10) 


y(t  =  0)  =  yo  (11) 


One  set  of  assumptions  that  has  been  hypothesized  (also 
proposed  by  H.  Weiss  [Ref.  26]  and  cited  by  Dolansky  [Ref.  9]) 
to  yield  equations  (10)  and  (11)  are  very  similar  to  the 
assumptions  listed  previously  for  the  square-law.   The  first 
three  assumptions  are  identical  to  the  first  three  square- 
law  assumptions. 

4(b).  Each  firer  knows  the  area  that  contains 
targets,  but  does  not  know  exact  target 
locations  or  the  consequences  of  his  fire. 

5(b).   Fire  from  surviving  units  is  distributed 
uniformly  over  the  target  area. 

The  state  solution  is  given  by 


b(xQ  -  x(t))  =  a(yQ  -  y(t))  (12) 
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hence  the  name  "linear-law"  attrition  process. 

The  time  solution  for  equations  (8)  and  (9)  is  given 
by  (see  G.  Weiss  [Ref.  35]) 


(yo-b/axQ) xQ  exp( (b/axQ-yo) at) 
X(t)  =  (yQ-b/axo  exp( (b/axo-yo) at)  (13) 


y(t)  =  yo-b/axQ  +  b/ax(t)  (14) 


Others  have  suggested  general  forms  of  homogeneous  force 
models  given  by 


|£  =  -  F(t,x,y)  (15) 


|£  =  -  G(t/x/y)  (16) 


For  example,  Willard  [Ref.  37]  considered  the  following 
Lanchester-type  equations: 


gf  =  "  ax°  f(t,x,y)  (17) 


|£  =  -  byC  f (t,x/Y)  (18) 

< 

Many  "general"  forms  contain  variable  coefficients  (e.g., 
a  =  a(t)  and  b  =  b(t)).   Taylor  [Ref.  30]  has  shown  that 
except  for  some  very  special  cases,  the  solutions  to  varia- 
ble coefficient  Lanchester-type  equations  are  very  complex 
and  in  many  cases  of  little  practical  use. 
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B.   STOCHASTIC  MODELS 

Koopman  [Ref.  22]  suggested  a  reformulation  of  the 
attrition  process  in  stochastic  form.   The  resulting  sto- 
chastic attrition  process  has  been  appropriately  called  by 
him  the  Lanchester  stochastic  process .   Others  have  subse- 
quently employed  a  stochastic  analysis  of  combat.   With  the 
inclusion  of  random  variations  in  the  attrition  process, 
analysts  have  attempted  to  better  represent  the  complex 
random  process  of  combat  with  the  goal  of  gaining  insights 
into  combat  not  available  from  deterministic  models. 

Snow  [Ref.  27]  shoAed  that  the  following  assumptions 
yield  the  square-law  attrition  process: 

(Al)   P(t+At,m,n:t,m+l,n)  =  anAt 

(A2)   P(t+At,m,n:t,m,n+1)  =  bmAt 

(A3)   Probability  of  more  than  one  casualty 
occurring  in  a  time  of  length  At  is  of 
the  order  of  magnitude  o(h)  where 

lim£(hL  =  o. 
h-0   h 

From  these,  he  developed  the  Chapman-Kolmogorov  forward 

equations : 

dP(t,m,n)  =  anp(t   +1   j  +  bmP(t,m,n+l)  -  (an+bm) P ( t,m, n) 
dt 

(19) 
If  one  assumes,  instead 

(Al)   P(t+At,m,n:t,m+l,n)  =  F(t,m+l,n)At 
(A2)   P(t+At,m,n:t,m,n+1)  =  G(t,m,n+l)At 
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(A3)   Probability  of  more  than  one  casualty 
occurring  in  a  time  of  length  At  is  of 
magnitude  o (h) . 


then  the  resulting  Chapman-Kolmogorov  forward  equations  are: 


dP(t^m,n)  =  F(t/m+1/n)p(tfin+1/n)  +  G(t,m,n+l)P(t,m,n+l) 
-  (F(t,m,n)  +  G(t,m,n))P(t,m,n)         (20) 


In  this  thesis,  this  process  will  be  called  the  general 
Lanchester-type  stochastic  attrition  process . 

Work  has  been  done  by  several  authors  to  obtain  general 
solutions  for  specific  stochastic  attrition  processes. 
Brown  [Ref.  7]  attempted  to  find  the  time  state  probabili- 
ties (P(t,m,n) )  for  a  class  of  homogeneous  force  stationary 
Markov  attritions  processes  (that  is,  when  F(t,m,n)  =  A(m,n) 
and  G(t,m,n)  =  B(m,n)).   He  concluded  that  the  expression 
derived  is  too  complicated  for  practical  use.   Isbell  and 
Marlow  [Ref.  14]  developed  a  general  solution  to  (20)  for 
a  square-law  stochastic  attrition  process.   If 
F(t,m,n)  =  an  +  3m  and  G(t,m,n)  =  bm  +  an  (Isbell  and  Marlow 
referred  to  a  and  3  as  operational  loss  rates) ,  with 
a  +  a  =  b  +  3 ,  then 


m  +n  -m-n 

P(t,m,n)  =  F(m,n:rn  ,n  )/ra  +n  \(exp  (  (b+3t) -1) 

o   o  |  o  o 

m+n  j  (21) 

exp(-(b+3) (m0+n0) t) 
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where  F(m,n:m  , nQ)  is  developed  from  a  recursive  relationship 


F  (111,11:111  ,n  )  = 


m  +an 

~'"^'  fuTnTz: — n — ; — \ *  (m,n:m  -±,n  ; 

o   o     (b+3)m  +  (a+ce)n         00 


n  +bm 

+  ,-,   7 Q, —, — ; — \ F(m,n:m  ,n  -1)   (22) 

(b+3)m  +(a+a)n      '    00 


where  F(m  ,n  :m  ,n  )  =  1. 
00   o   o 

Clark  [Ref.  8]  developed  a  general  solution  for  the  time 
state  probabilities  for  the  linear-law  stochastic  attrition 
process  (i.e.,  F(t,m,n)  = amn  and  G(t,m,n)  =  bmn) ; 


m 


n   ! ,  , . k-n+j-m   o 
o   o   (-1)         a 

P(t,m,n)  =   Z    I 


m_-m  n  -n 

b  °    (m  )  !  (n  )  ! 
o     o 


.    .  n  -n+m  -m 

(a+b)  °    °   mi  n!  (k-n)  !  ( j-m)  I  (n  -k)  ! 


m  -j 


exp (- (a+b) jkt 


\ 


(23) 

(However,  he  did  not  compare  his  results  with  those  for  the 
probability  of  winning  (i.e.,  "true"  absorption  probabili- 
ties) .   This  expression  is  too  complicated  to  provide  any 
insights  by  itself.   This  author  has  no  knowledge  of  any 
other  general  solutions  developed  for  the  time  state  proba- 
bilities.  It  is,  however,  a  straight  forward  task  to  compute 

P(t,m  ,n)  and  P(t,m,n  ),  which  are  used  in  the  included 
'  o  o 

computer  program  as  a  partial  check  for  accuracy. 
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Work  has  also  been  done  to  develop  exact  expressions  for 
the  probability  of  winning,  for  a  stationary  Markov  attrition 
process.   The  exact  expression  for  the  probability  of  winning 
for  the  linear-law  stochastic  attrition  process  is  easily 
obtained  from  classical  random  walk  results  [Ref.  28]: 

n  n    M       ,       m       n  -n 

o  o   m  +n  -n-1        o   .    o 

PY  =  l      P(o,n:mo,no)  =  I   (  °  ?    )  ,-^)    <_*-> 

n=l  n=l     o 

(24) 

Smith  [Ref.  26]  extended  work  done  by  Brown  [Ref.  7]  and 
developed  the  explicit  expression  for  the  square-law  stochas- 
tic attrition  process: 

m 
o 

Pv  =   S   P(m,0:rn  ,n  )  (25) 

X    m=l         °   ° 

where 

.     m      m  -j   m  +n  -m-1 

bm     o  ,  , »  o  J  .o   o     twi_/  ■  ,  t  \ 

P(m  0-m  n  )  =  T 2-1  7    r("1)      3 r  (b/a  y+1) 

^m,u.mo,no}    Lan  +bm  J  I    [  _  a  j+1)(j_m+1)J 

o    o   j =m    oJ        o      J     J 

(26) 

This  is  an  impressive  looking  expression,  but  from  a  practical 

viewpoint,  it  is  not  very  useful.   A  computer  is  needed  to 

solve  this  equation  for  any  appreciable  force  sizes,  but 

one  immediately  runs  into  overflow  and  underflow  problems. 

To  solve  these  complications,  recursive  relationships  must 

be  used;  but  if  one  looks  carefully  at  the  relationships  . 

used,  he  is  lead  back  to  the  initial  recursive  relationships 

used  to  obtain  the  solution: 
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A(mQ,no) 

P(0,n:m  ,n  )  =  —, c  ,  _  , r  P(0,n:m  -l,n  ) 

o   o    A(m  ,n  )+B(m  ,n  )     '    o    o 
o  o     o  o 

B(m  ,n  ) 

+  T7 TT5T r  P(0,n:m  ,n  -1)    (27) 

A(mo/no) +B (m  ,nQ)         o'  o 


It  was  found  by  this  author  that  the  most  efficient  method 
to  solve  the  stochastic  model  for  the  probability  of  winning 
was  by  using  the  above  recursive  relationship,  or  one  devel- 
oped by  Springall  [Ref.  28],  not  the  general  force  level 
solutions  (See  Appendix  A) . 

Springall  developed  a  recursive  relationship  .for  the 
distribution  of  the  time  of  battle  termination: 

n 
.r,  _   °    .,_  .,  _%„(r) 


E(t  )  =   I      A(%  +l.n)3C  '+1 
n=n,  +1    ^       dp 


*bp  '  -  ^p 

m 

°  (r) 

+   z     B(ra,a  +l)Xj>ri   J..  (28) 


m=rab 


B(m,n,  +1)XV  '   .- 
+  1     ^P    m,n^  +1 

P 


where 


x(0)   = i (29) 

m  ,n    A(m  ,n  )+B(m  ,  n  ) 
o   o      o '  o     o     o 

Blm.n+Dxf1., 

x(0)    =  o  (30) 

m  ,n    A(m  ,n)+B(m  ,n) 

A(m+l,n  )X^\ 
/nv  o   m+l,n 

x(°)    -  o  (3D 

m,no    A(m,n  ) +B (m,n  ) 
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x(0)  m+l,n  m,n+l        ._,,> 

m,n  A(m,n)+B(m,n)  K      ' 


for  i  >  1 


ixti-D 

/•i  m  , n 

m  ,n    A(m  ,n  ) +B (m  ,n  )  v   ' 

o   o      o  o     o  o 


ix(i-l)+B(m  n+1)xd) 
/•x       in  ,n     o      m  ,n+l 

mQ,n    A(mo/n)+B(mo,n)  K°*' 


iX(l"1)+A(m+l/n^)X(iJ 
/ • v        m,n          o   m+l,n 
Y(i)    =       o     __o  ,-,-> 

m,n     A(m,n  )+B(m,n  )  v   ; 

o         o       o 


Y(i)      iXf[1n1)+A(m+lfn)Xi!i^   +B  (m,n+l)  X*1* 
X      =   m , n  m+ 1 , n  m, n+1 


m,n 


A(m,n)  +B  (m,n) 

(36) 


From  studying  the  equations  presented  in  this  section,  one 
gains  an  appreciation  for  the  increase  in  complexity  of  the 
model  by  the  inclusion  of  uncertainty  in  the  attrition 
process.   It  is  apparent  that,  although  exact  solutions  may 
be  obtained  for  some  stochastic  models,  numerical  results 
are  not  readily  generated  by  hand.   A  computer  is  therefore 
essential  for  the  numerical  solution  of  the  simple  stochastic 
models . 

It  is  surprising  that  more  use  has  not  been  made  of  the 
computer  in  investigating  stochastic  models,  especially  the 
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use  of  computer  graphics.   For  example,  the  five  following 
figures  are  plots  of  the  time  state  probabilities  for  a 
stochastic  model  at  3,  25,  50,  75,  and  100%  of  the  time  to 
completion  of  the  equivalent  deterministic  model.   It  is 
believed  that  something  of  this  nature  has  not  been  utilized 
in  the  open  literature  before.   It  gives  some  interesting 
insights  into  the  progress  of  a  stochastic  model.   The  state 
probability  seems  to  start  as  a  spike  at  (m  , n  )  and  pro- 
gress{' over  time  somewhat  as  a  drop  of  water  would.   As  time 
increases,  the  drop  spreads  and  moves  toward  the  banks 
(boundaries  of  the  state  space) .   The  angle  at  which  the 
wave  moves  is  determined  by  the  relationship  between  the 
attrition  rates  of  the  process  (in  this  case,  the  model  is 
the  square-law  stochastic  attrition  process,  so  the  rela- 
tionship is  between  a,  b,  m,  and  n) .   As  the  "wave"  hits 
the  wall,  it  "sticks"  (is  absorbed)  ,  thus  giving,  as  t-*-00, 
the  distribution  of  win  probabilities.   (It  is  interesting 
to  note  that  this  relates  the  work  of  Clark  [Ref.  8], 
Weiss  [Ref.  35],  Morse  and  Kimball  [Ref.  22],  and  Snow 
[Ref.  27] ,  who  were  concerned  with  the  solution  to  the 
Chapman-Kolmogorov,  to  the  work  of  Smith  [Ref.  26],  Brown 
[Ref.  7],  etc.,  who  were  concerned  with  the  win  probabilities) 
Note  that  in  the  present  case,  Y  has  a  very  large  probability 
of  winning.   If  the  forces  were  moved  closer  to  parity,  more 
of  the  wave  would  be  absorbed  by  the  Y  force  boundary, 
giving  X  a  larger  probability  of  winning. 
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Plot  of  P(t,m,n)  for  Fixed  t 
FIGURE  3 
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The  above  discussion  has  given  a  description  of  what 
occurs  with  the  probability  mass  over  time.   Another  des- 
cription .of  the  process  is,  it  acts  as  a  combination  of 
diffusion  and  convective  transport  of  probability  mass. 
This  concept  enhances  ones  understanding  of  the  nature  of 
so-called  diffusion  approximations  of  Markov  processes 
(i.e.,  state  space  assumed  continuous).   Additionally,  it 
provides  insights  into  temporal  changes  in  variance  in  force 
levels  (see  Chapter  III,  Section  C) .   Thus,  it  is  felt  that 
computer  graphics,  if  utilized  more  fully  could  provide 
easily  obtainable  insights  into  the  dynamics  of  a  combat 
model,  that  would  be  difficult  to  obtain  otherwise. 

C.   COMPARISONS 

Two  types  of  comparisons  have  been  made  between  stochas- 
tic and  deterministic  models;  probability  of  winning,  and 
the  expected  force  levels  against  the  time  history  of  the 
force  levels  of  the  deterministic  model.   Two  recent  studies 
have  been  done  comparing  the  probability  of  winning  with 
the  certainty  of  winning  or  losing  in  the  deterministic 
model  [Ref.  20  and  Ref .  28].   Both  came  up  with  basically 
the  same  conclusions. 

(1)  As  force  levels  increase,  the  probability 

of  winning  approaches  the  certain  probability 
of  winning  or  losing  in  the  deterministic 
model. 

(2)  Significant  differences  occur  only  when  the 
forces  are  near  parity. 
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Lee  and  Wannasilpa  [Ref .  20]  came  to  their  conclusions  for 
the  square-law,  linear-law,  and  mixed-law.   Springall  arrived 
at  the  same  conclusions  for  a  slightly  more  complicated  model, 

Snow  [Ref.  27]  showed  that  the  time  history  of  the 
expected  force  levels  differed  from  the  time  history  of 
the  deterministic  force  levels;  that  is  m(t)-x(t)  ^   0. 
Clark  [Ref.  8]  defined  this  as  bias .   This  seems  to  be  an 
appropriate  term  and  will  be  used  here.   Snow  apparently 
first  developed  an  explicit  expression  for  this  bias  for 
the  square- law; 

n 

dm      —       ° 

^r  =  -  an  +  a  I   P(t,0,n)  (37) 

n=l 

Clark  [Ref.  8]  empirically  showed  that  bias  exists  for  the 
linear-law.   All  researchers  seem  to  have  arrived  at  the 
same  general  conclusions: 

(1)  As  initial  force  levels  increase,  the  bias 
decreases. 

(2)  The  longer  a  battle  lasts,  the  greater 
the  bias. 

Powers  and  Taylor  [Ref.  24]  and  Hanna  [Ref.  12]  have 
compared  differences  in  optimal  time  sequential  fire  distri- 
bution policies  for  deterministic  and  stochastic  models: 
Hanna' s  investigation  indicates  that  there  is  a  real  differ- 
ence between  the  two  (at  least  for  small  numbers  of  total 
combatants)  . 
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Thus,  some  work  has  been  done  to  investigate  the  differ- 
ences between  deterministic  and  stochastic  models,  especially 
in  the  areas  of  who  is  going  to  win  and  the  average  time 
histories  of  the  force  levels.   However,  the  analysis  done 
up  to  the  present  has  been  somewhat  limited  in  scope  and 
does  not  allow  an  analyst  to  say  with  any  degree  of  cer- 
tainty, when  the  complex  random  process  of  combat  can  be 
adequately  represented  by  a  deterministic  model. 
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III.   COMPARISON  OF  DETERMINISTIC  AND  STOCHASTIC 
MODELS  FOR  "SQUARE  LAW"  ATTRITION  PROCESS 

It  would  be  most  useful  to  be  able  to  make  comparisons 
between  stochastic  and  deterministic  models  used  today. 
This  is  not  feasible  within  the  scope  of  this  paper  for 
several  reasons. 

1.  It  is  not  clear,  in  many  cases,  when  the  two  models 
are  equivalent.   For  example,  one  could  try  and 
compare  DYNTACS-X  and  Bonder/I .U. A. ,  but  before 
this  comparison  could  be  made,  one  would  have  to 
insure  that  the  input  parameters  are  equivalent. 

In  the  case  of  these  two  models,  that  is  not  a 
trivial  task,  as  the  structure  of  the  two  models 
is  quite  different,  particularly  in  route  selection. 
About  the  best  that  can  be  done  is  to  construct  the 
input  parameters  for  each  model  and  assume  they 
are  equivalent. 

2.  Parametric  analysis  is  not  feasible  in  a  complicated 
model.   Complicated  models,  both  stochastic  and 
deterministic,  are  expensive  to  run.   Because  of 
this  and  the  complex  interactions  among  the  variables 
of  a  model,  parametric  analysis  is  not  only 
prohibitively  expensive,  but  difficult  to  do. 

3.  It  is  difficult  to  gain  insights  in  a  comparison 
of  simple  stochastic  and  deterministic  models,  and 
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essentially  impossible  to  gain  many  useful  insights 
in  a  comparison  of  complicated  models. 
Lanchester-type  models  may  be  thought  of  as  resulting  from 
an  aggregation  of  the  myriad  of  combat  environmental  varia- 
bles into  constant  attrition  rate  coefficients  and  time 
varying  force  levels  and,  moreover,  as  complementary  complex 
system  models.   Such  idealizations  will  be  considered  in 
this  thesis  with  the  knowledge  they  are  idealizations,  but 
that  they  may  provide  insights  that  can  be  used  as  points 
of  departure  for  the  investigation  of  more  complex  models 
of  combat.   Specifically,  two  models  will  be  compared;  the 
deterministic  Lanchester  square-law  and  the  square-law 
stochastic  attrition  process.   It  is  believed  that  trends 
found  from  this  investigation  may  be  carefully  used  to 
determine  trends  for  a  comparison  of  more  complex  models. 
As  was  shown  in  the  literature  survey,  the  square- law 
is  a  special  case  of  the  general  differential  equations 


jnr  -  -  F(t,x,y)  (15) 


dt 


g£  =  -  G(t,x,y)  (16) 


where 


F(t,x,y)  =  ay 


and 
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G(t,x,y)  =  bx 


giving 


dx  ,  -, , 

at  =  "  ay  (1> 


£  =  -  bx  (2) 


with  initial  conditions 


X(t  =0)  =  xQ  (3) 


Y(t  =  0)  =  yQ  (4) 


Also  shown  was  the  general  form  of  the  stochastic  model: 


dP(tdt'n)  =  F(t'm+1'n)p(t'm+1'n>  +  G(t,m,n+l)P(t,m,n+l) 


-  (F(t,m,n)  +  G(t,m,n) )P(t,m,n)       (20) 


with  boundary  conditions 


dP  ( t ,  m.  ,  n   ) 

d°P  bP  =  F(t.^,+l.nb  )P(t,^p+l,nb^) 


+  G(t'mbp'nbP+1)P(t'mbp'nbp+1)      <38) 
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dP(t,m,n,  ) 

dt   P?   =  F(t,m+l,nbp)P(t,m+l,n   ) 


+  G(t,m,nb  +l)P(t,m,nb  +1) 


F(t,m/nb  )P(t,m,nb  ) 


m  >  m 


np 


(39) 


dP(t,m   ,n) 

■^2 =  F(t,m  +l,n)P(t,m  +l,n) 

np  np 


dt 


+  G(t,m   ,n+l)P(t,m   ,n+l) 
np  np 


-  G(t,m   ,n)P(t,m   ,n) 
np        np 

n  >  "bp 


(40) 


and  initial  conditions 


P(0,m,n) 


1   m  =  m   and  n  =  n 
o  o 

0   otherwise 


(41) 


Substituting  the  square-law  attrition  rates  gives  the 
square-law  stochastic  attrition  process: 


dP  (t,m,n) 
dt 


anP(t,m+l,n)  +  bmP ( t,m,n+l)  -  (an+bm) P ( t,m,n) 

(19) 


and  boundary  conditions 
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dP  ( t ,  m   ,  n,  ) 

dT* =   anbpP(t'ntp+1'nbP)  +  bVP(t/mbp'V+1) 

(41) 


dP(t,m,n.   ) 

dt 


=  ar^  P(t,m+l,n   )  +  bmP(t,mfnb  +1) 
-  anbpP(t/m/nbp) 


(42) 


dP(t,m   ,n) 
dt 


=  anP(t,m  +l,n)  +  bir^  PU,!^  ,n+l) 


bm,  P(t,m   ,11) 
dp     np 


(43) 


and  initial  conditions 


P(o,m,n) 


( 


m  =  m   and   n  =  n 
o  o 

0   otherwise 


(44) 


It  was  noted  that  the  square-law  comparisons  have  been 
made  previously.   However,  as  far  as  the  author  can  tell, 
no  previous  comparisons  have  been  made  contrasting  the  models 
for  breakpoint  rules  other  than  annihilation.   Taylor  [Ref .  30] 
has  pointed  out  (and  Adkins  [Ref.  1]  examined  in  further 
detail)  that  a  neglected  area  in  the  modeling  of  land  combat 
is  the  modeling  of  battle  termination.   In  actual  battles, 
several  different  events  might  cause  battle  termination. 
If  one  of  the  opposing  forces  is  annihilated,  then  the  battle 
ends;  however,  this  is  a  rare  event  essentially  never  observed 
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in  real  combat,  with  a  few  exceptions  (e.g.,  Iwo  Jima,  Alamo, 
etc.).   If  one  of  the  opposing  forces  surrenders  unconditionally, 
the  battle  will  also  terminate.   The  third  and  most  common 
event  which  will  result  in  battle  termination  is  that  of  one 
of  the  opposing  force  breaking  contact  with  the  enemy  and 
withdrawing  from  the  battlefield  [Ref.  13].   There  are  other 
possible  events  which  might  result  in  battle  termination, 
but  they  will  not  be  considered. 

For  the  purposes  of  further  discussion,  a  breakpoint  is 
defined  to  be  that  state  of  a  battle  at  which  a  unit  con- 
siders itself  no  longer  capable  of  performing  its  mission 
and  as  a  result  elects  to  break  contact  with  the  enemy  and 
withdraw  from  the  battlefield.   Therefore,  when  a  unit  with- 
draws from  the  battlefield  strictly  to  avoid  further  combat, 
the  unit  is  considered  to  have  reached  its  breakpoint  and 
lost  the  engagement. 

In  many  current  land  combat  models  a  unit's  breakpoint 
is  determined  by  the  percentage  of  casualties  sustained  by 
that  unit  [Ref.  1] .   Although  available  empirical  evidence 
shows  that  a  deterministic  battle  termination  model  does  not 
appear  reasonable,  it  is  an  adequate  "first  cut"  [Ref.  13]. 
And  a  breakpoint  determined  by  the  percent  of  casualties 
sustained  is  more  reasonable  than  one  that  assumes  the  battle 
goes  to  annihilation.   For  this  reason,  the  investigation 
was  done  with  deterministic  breakpoints  determined  by  the 
percent  of  casualties. 
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A.   PROBABILITY  OF  WINNING 

Hypothesis  1-1 

i 
As  the  breakpoint  force  levels,  m,   and  n,   are 

moved  closer  to  the  initial  force  levels,  m   and 

o 
n  ,  respectively;  the  difference  m  the 

probability  of  winning  between  the  two  models 

increases . 

Lee  and  Wannasilpa  [Ref.  20]  deduced  that  as  the  initial 
force  levels  are  increased,  the  difference  between  the  proba- 
bility of  winning  in  the  deterministic  model  and  the  stochas- 
tic model  decreases.   In  fact  it  appeared  that  the  difference 
disappears  in  the  limit.   However,  the  conclusion  reached 
by  Lee  and  Wannasilpa  was  based  only  on  battles  to  annihila- 
tion.  A  look  at  the  recursive  relationship  used  in  Appendix 
A  to  solve  for  the  probability  of  winning  is  based  on  a 
summation  from  m,  +1  to  m  .   Furthermore,  the  recursive 
relationships  are  from  itl   to  m  .   These  seem  to  indicate 
that  the  differences  in  the  probability  of  winning  decrease 
not  just  as  m  increases,  but  as  m  -m,   increases.   The 


o 


o  ™bp 


following  figure  seems  to  support  this.   Note,  that  for  a 
deterministic  model,  the  probability  that  Y  wins,  Py ,  has 

a  discontinuous  jump  from  one  to  zero  as  the  ratio 

bmo2(l-(l-Fx)2) 

« 7T—   passes  through  one.   (Note:   Fv  and  F   are 

ano2(l-(l-FY)2)  X      Y 

the  fractions  of  initial  combatants  each 

side  is  willing  to  lose.)   For  the  square- law  stochastic 

attrition  process,  however,  P   moves  continuously  from  one 

to  zero  reaching  about  0.5  at  a  ratio  of  one.   Note  that  as 

Fv  and  F..  decrease,  m,   and  n,   increase,  so  m  -m,   and 
X      Y  bp      bp  o  bp 
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Py   for    Varying   Breakpoints 


a=0.008 
b=0.004 
No=10 


n  -n,   decrease.   As  they  decrease.  P.,  for  the  stochastic 
o  bp  *  '      Y 

model  moves  further  and  further  from  P   for  the  deterministic 
model.   This  was  found  to  be  true  for  all  combinations  of 
parameters  considered. 

Hypothesis  1-2 

For  fixed  breakpoint  casualty/initial  force 
ratios  (F   and  F  ) ,  the  difference  between 
the  probability  of  winning  for  the  deter- 
ministic model  and  the  stochastic  model 
decreases  as  the  initial  force  levels  increase. 

Lee  and  Wannasilpa  [Ref.  20]  deduced  the  above  except 

their  conclusion  was  restricted  to  battles  to  annihilation. 

The  above  hypothesis  generalizes  the  conclusion  to  include 

any  set  of  fixed  F   and  Fv .   Once  hypothesis  1-1  along  with 

A  X 

Lee  and  Wannasilpa1 s  conclusion  are  accepted,  hypothesis 
1-2  follows.   The  following  figures  support  this.   It  is 
interesting  to  note  that  the  plots  of  Py  are  similar  for 
all  fixed  sets  of  F   and  F  .   The  only  difference  is  in  the 
angles  of  the  curves.   As  F„  and  FY  increase,  the  result 
is  similar  to  rotating  the  curves  counter-clockwise.   This 
causes  an  increase  in  the  difference  between  P„  for  the  two 
models . 

Hypothesis  1-3 

As  the  forces  move  away  from  parity,  the 

difference  between  the  probability  of 

winning  for  the  two  models  becomes  negligible. 

This  was  the  final  conclusion  reached  by  Lee  and  Wanna- 
silpa.  Although  their  conclusion  was  for  battles  to 
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Py    for   Fixed    Fx  and    Fy 


a=0.008 
b=0.00^ 

Fxf1-0 

Fy=1.0 


bM„2(l-(l-FY)2) 


figl\!;:  7 
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'v  for  Fixed  Px  and  Py 


1.0 


P(Y  wins) 


a=0.008 
b=0.004 
Fx=0.8 

Fy=0.8 


aN0*91-(l-Fy)*) 
FIGURE  8 
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l.g(Y  wins) 


Py   for    Fixed    F^  and    F„ 


a=0.008 
b=0.004 
Fx=0.6 

Fy=0.6 


aNo^Cl-Cl-Kyr) 

fig  lira-:  9 
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PY   for    Fixed    P„  and    Py 


0.2- 


0.1- 


0.0 


a=0.008 
b=0.00^ 
Fx=0.4 

Fy=0.4 


bM^(l~(l-I-V)-) 


FJOllRE    10 
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P      for   Fixed    Fx  and    Fy 


P(Y  vd.ns) 
1.0- 


0.0-- 


bH^Cl-q-r^) 
aN0^(l-(l-Fv)^) 


FI GU  RE    11 
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annihilation,  the  present  findings  do  not  alter  the  conclu- 
sion for  other  fixed  force  level  breakpoints.   The  only 
change  is  in  how  far  from  parity  the  forces  must  be  to  make 
the  difference  negligible. 

The  general  consequence  of  this  hypothesis  is  that  unless 
the  forces  are  near  parity,  there  is  little  difference  in 
the  probability  of  winning  between  the  two  models.   From 
this  viewpoint,  then,  if  one  side  is  going  to  win  decisively, 
there  is  little  reason  to  use  the  more  expensive  stochastic 
model.   Indeed,  this  has  some  intuitive  appeal:   The  more 
decisively  a  side  wins  the  battle,  the  less  influence  any 
uncertainties  will  have  on  the  outcome.   On  the  other  hand, 
if  the  forces  are  near  parity,  the  uncertainties  in  the 
evolution  of  the  dynamics  of  combat  may  cause  the  expected 
loser  to  win. 

If  the  previous  hypotheses  are  combined  with  this  hypothe- 
sis, some  interesting  trends  can  be  observed.   It  has  been 

concluded  that,  with  increasing  m  -m,   and  n  -n,  ,  the 

3   o  bp      o  bp 

differences  in  the  probability  of  winning  between  the  two 
models  decreases.   When  combined  with  hypothesis  1-3,  this 
says  that,  with  casualties  on  each  side  on  the  order  of  20 
or  more,  and  with  forces  not  near  parity,  there  is  no 
significant  difference  between  the  probability  of  winning 
for  each  model.   Even  if  the  forces  are  near  parity,  if  the 
force  levels  are  very  large,  and  both  sides  are  willing  to 
take  20  or  more  casualties,  the  difference  in  model  win 
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probabilities  is  negligible.   In  these  cases,  if  only  the 
win  probabilities  were  considered  as  measurement  criteria 
between  the  two  models,  one  would  use  the  less  expensive 
deterministic  model  as  an  acceptable  substitute  for  the 
more  realistic  stochastic  model. 

One  must  be  cautious  of  this  conclusion.   It  would  be 
inaccurate  to  make  this  conclusion  unless  the  only  objective 
of  the  model  is  to  determine  which  side  will  win.   In  most 
cases,  there  are  other  insights  analysts  try  to  gain  from 
the  model,  such  as  how  many  casualties  does  each  side  take, 
at  what  points  in  the  battle  are  the  most  casualties  being 
taken,  how  long  will  the  battle  last,  etc. 

Nevertheless,  the  conclusions  of  this  section  do  give 

indications  of  when  the  two  models  will  not  give  the  same 

winners.   In  these  cases,  no  further  evaluation  of  the  models 

is  needed,  because  it  is  obvious  the  deterministic  model  is 

not  an  acceptable  substitute.   As  Springall  [Ref.  28]  said; 

Probably  the  single  most  important 
comparison  is  between  the  two  predic- 
tions of  which  side  is  going  to  be 
victorious.   If  the  two  methods  cannot 
agree  on  this,  there  is  little  hope 
for  agreement  on  the  subsidiary 
attributes . 

B.   TIME  HISTORY  OF  EXPECTED  FORCE  LEVELS 

If  a  determination  has  been  made  that  the  models  will 
predict  the  same  winner,  further  analysis  may  be  made  to 
see  how  well  the  models  compare.   Clark  [Ref.  8]  suggests 
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that  a  useful  comparison  is  the  time  history  of  the  expected 

number  of  survivors.   In  many  models,  one  of  the  most  utilized 

variables  (performance  and/or  proxy  variables)  in  measures 

of  effectiveness  for  systems  being  evaluated  by  the  models 

is  the  number  of  survivors.   It  has  been  stated  previously 

that  a  bias  exists  in  the  square-law  and  the  linear-law. 

Lanchester  himself  suspected  as  much;  he  recognized  that 

his  differential  equations  were  approximations  to  the  casualty 

rates  likely  to  be  experienced  in  an  actual  battle.   Quoting 

from  Lanchester  [Ref.  19]: 

Since  the  forces  actually  consist  of  a 
finite  number  of  finite  units  (instead  of 
an  infinite  number  of  infinitesimal  units) 
the  end  of  the  curve  must  show  discontinuity, 
and  break  off  abruptly  when  the  last  man  is 
reached;  the  law  based  on  averages  evidently 
does  not  hold  rigidly  when  the  numbers  become 
small. 

Lanchester  stated  that  his  differential  equations  are  based 

on  averages  implying  an  underlying  stochastic  process.   He 

also  suggests  that  his  differential  equations  may  be  good 

approximations  only  so  long  as  the  force  sizes  are  large. 

Although  his  reasoning  is  heuristic,  Snow  [Ref.  27]  verified 

that  the  square-law  is  biased  for  small  numbers.   He  showed 

that 


n 

-,—     _      o 

~  =  -an  +  a   E   nP(t,0,n)  (37) 

n=l 


m 

dn      —       ° 

~  =  -bm  +  b   Z   mP(t,m,o)  (45) 

m=l 
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The  conclusion  is  that  the  solutions  to  Lanchester's  equa- 
tions for  the  square- lav;  are  good  approximations  to  expected 
survivors  so  long  as 

n 
o 

a   I   nP(t,0,n)  2  0  (46) 

n=l 


m 
o 

b   E   mP(t,m,0)  =  0  (47) 

m=l 


Intuitively  this  seems  sound.   (46)  and  (47)  express  the 
boundary  effects  at  time  t.   At  early  times  in  the  battle, 
the  expected  force  levels  are  near  the  initial  force  levels, 
and  there  is  very  little  of  the  probability  mass  at  the 
boundary.   Thus  the  boundary  effects  are  small.   As  the 
battle  nears  termination,  though,  much  of  the  probability 
mass  is  near  the  boundary  and  (46)  and  (47)  will  be  signi- 
ficant enough  to  effect  the  process.   As  there  are  no  boun- 
dary effects  in  the  deterministic  process,  there  will  be 
a  difference  in  force  level  trajectories  of  the  two  processes. 

Clark  provided  a  derivation  of  (37)  and  (45).   Additionally, 
he  provided  the  following  argument  showing  that  the  linear- 
law,  which  at  first  appears  to  be  unbiased,  is  not.   Through 
a  simple  derivation,  it  can  be  shown  for  the  linear-lav/: 


^  =  -  aE(MN:t)  (48) 
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Il  =  -  bE(MN:t)  (49) 


There  are  no  bias  terms  in  this  expression  (as  opposed  to 
(37)  and  (45))  since  the  random  variable  MN  has  a  value  of 
zero  whenever  M  =  0  or  N  =  0 .   Consequently,  the  assumed 
differential  equations  for  the  linear-law  appear  unbiased. 
However,  even  though  the  analysis  of  the  solutions  to 
these  equations  show  that  the  state  solution  is  unbiased, 
the  time  solution  interjects  a  source  of  bias.   This  is 
easily  seen  as  a  result  of  elementary  probability  theory: 


E(MN)  =  E(M)-E(N)   if  an  only  if  M  and  N 

are  independent. 


In  the  linear-law,  they  are  obviously  dependent  on  each 
other.   The  rate  of  change  of  each  force  level  is  dependent 
on  the  product  of  the  two  force  levels,  causing  the  force 
levels  to  be  dependent. 

Clark  [Ref.  8]  also  gave  numerical  examples  to  show  the 
amount  of  bias  resulting  from  the  difference  in  models.   His 
results  showed  that  there  are  significant  amounts  of  bias 
resulting  in  both  the  linear-law  and  square- law  cases. 
There  are  two  things  that  should  be  noted  about  his  results: 
First,  his  results  are  for  small  numbers  of  combatants;  and 
second,  his  results  are  only  for  battles  to  annihilation. 
It  was  not  shown  in  his  dissertation  that  there  is  a 
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significant  bias  with  large  numbers  of  combatants  or  with 
breakpoints  other  than  annihilation. 

The  following  hypotheses  are  based  on  extension  of  his 
conclusions  to  include  larger  forces  and  fixed  force  level 
breakpoints.   Appendix  D  is  a  derivation  of  the  square-law 
stochastic  attrition  process  showing  that  the  bias  noted 
by  Snow  exists  with  any  set  of  fixed  breakpoint  force  levels 
Unfortunately,  the  form  of  the  bias  is  not  as  clear  as  that 
for  annihilation,  since  there  are  two  additional  terms.   The 
use  of  the  graphic  aids  noted  earlier  gives  some  indication 
that  non-zero  breakpoint  force  levels  do  not  change  the 
general  progress  of  the  battle.   A  look  at  the  following 
figures  will  show  that  the  time  state  probabilities  progress 
in  the  same  manner  as  in  the  case  of  battles  to  annihilation 
Figures  12-15  are  plots  of  the  time  state  probabilities  of 
the  same  model  shown  earlier  except  the  breakpoint  force 
levels  are  eight  instead  of  zero  (annihilation) .   Figures 
16-19  are  the  same  except  the  breakpoint  force  levels  are 
sixteen.   The  differences  are  that  the  "wave"  is  absorbed 
by  the  boundary  at  (m,  ,n)  and  (m,n,  )  instead  of  at  (0,n) 
and  (m,0) .   The  question  now  is,  What  is  the  effect  on 
the  bias? 

A  first  analysis  of  (101)  and  (102)  might  give  the 
indication  that  the  bias  is  always  positive  (i.e., 
m(t)-x(t)>  0  for  all  t  >  0) .   However,  this  is  not  true, 
as  can  be  seen  in  the  following  argument. 
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a=0.OO8 
b=0.00^ 
Mo=<*0 
Mbp=8 

"bt>=8 


Survivors 


Survivors 


Plot   of    P(t,m,n)    for    Fixed    t 
FIGURE     12 
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P(t,m,n) 


•=0.008 

b=0.004 

Mo=40 

Mbp=8 

No=40 

"bp=8 

t=60 


Survivors 


Plot   of   P(t,m,n)    for   Fixed    t 
FIGURE   13 
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P(t,m,n) 


a=0.008 
b=0.004 
Mo=40 

mEho 

Nbn=8 


Survivors 


Plot   of   P(t,m,n)    for   Fixed    t 
FIGURE     14 
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P(t,m,n) 


•=0.008 

b=0.004 

Mo=40 

Mbp=8 

N0=V) 

Nbp=8 
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Define 


Ax   =   m(t)    -    x(t)  (50) 


Ay   =   n(t)    -   y(t)  (51) 


then 


Av(t=0)    =   m(t=0)-x(t=0)    =   -    an      +    an      =0  (52) 

2\  O  O 


Ay(t=0)    =   n(t=0)-y(t=0)    =   -   bmQ   +  bmQ   =    0  (53) 


Let 


n 

S.r(t)    =   a[    Jp    nP(t,m      ,n)    +   n,     m.     P(t,m      +  l,n,     ) 
Y  _  '    np  Dp  Dp         '    np  bp 

n  nbp 

m 

+  v/m  +1p(t'm'V]  (54) 

c   m=m      +1  c 

np 

m 

S    (t)    =  b[    S      nfl?<t,m,iL     )    +  i^  n      P(t,m      ,n      +1) 
m=m  ^  tr       tr  l 

np 

n 
o 

+  ni  E  P(t,m.      ,n)  ]  (55) 

np  ,  ,  op 

^  n=n,     +1  c 

bp 


If 
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dx  , .. . 

at  ■  -  «y  (1) 


and 


^  =  -  an  +  Sy(t)  (101) 


then 


dm   dx      —     ... 
__-__=_  an  +  sY(t)  +  ay 


d(m-x)  _     ,-  .      ,  . 
— -^ a(n-y)  +  SY(t) 


Substituting  (50) 


dAy 

-^  =  -  aAy  +  Sy(t)  (56) 


and  equivalently 


dA 

-gi-  -  bAx  +  SY(t)  (57) 


with  initial  conditions 


Ax(t  =  0)  =  0  (52) 
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AY(t  =  0)  =  0 


(53) 


Define 


A  = 


t*X  \ 


and   U  = 


8,(1) 
lSX(t> 


then 


di 
dt 


0   a 


A+  U(t) 


lb   o 


(58) 


and  the  solution  of  (58)  is  given  by 


where 


A(t)  =  e  AtA(t=0)  +   /  e"A(t"T)U(x)  dT   (59) 


A  = 


0   a^ 
b  o\ 


Proof: 


dt 


■Ae    A(t=0)  +  U(t) 


A   /  e"A(t~T)U(T)dT  -  AA(t) 


-AA(t)  +  U(t) 
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from  Sylvester's   Theorem 


|  cosh    ( "\  ab    t)  -\la/b   sinh(^v  ab    t) 

e  = 


-"\,  b/a  sinh(^  ab    t)  cosh    (~\j  afo   t) 


Hence 


t 
/ 
0 


Ax(t)    =      /   S    (x)cosh(\'ab    (t-x))    -   Sx(t)^  a/b    sinhtyab    (t-x))dx 

(60) 


Ay(t)  =   /  Sx(x)cosh(  YaF  (t-x)  )  -  S^Tj^b/a  sinh  Cyab~  ( t-x)  )  di 
0  (61) 


If  one  now  takes  the  case  where  one  force  is  annihilated 

(m   =  n,   =  0)  .  then  Sv(t)  and  S.7(t)  simplify  to 
np    bp      '        X         Y       c        JL 

n 
o 

S„(t)  =  I      nP(t,0,n)  (62) 

*       n=0 

m 
o 

Sv(t)  =   Z   mP(t,m,0)  (63) 

X      m=0 


which  are  simply  the  average  number  of  Y  survivors  when  X 
is  annihilated  at  time  t  and  the  average  number  of  X  sur- 
vivors when  Y  is  annihilated  at  time  t  respectively.   It 


Sylvester's  Theorem:   For  a  polynomial  function  or 
convergent  power  series  of  a  matrix,  when  the  matrix  has 
distinct  eigenvalues,  A., 

n        n3  A-A.I 

F(A)  =   Z  F(A.)   Z  (-r J — )  where  A  is  an  nxn  matrix 

i-i     x    j=i  VAj 
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seems  intuitively  appealing  that  if  X  wins  decisively, 

Sv(t)  -*■  0 .   If  this  is  true,  then  (60)  will  become  negative 

and  (61)  will  become  positive.   This  seems  to  be  confirmed 

by  numerical  results  obtained. 

Another  interesting  case  to  examine  is  the  case  of 

complete  parity  (a=b,  m  =n  ,  m  =n,  )  where  m  =n,  =0. 
r     *  2  '      o   o'   np  bp         np  op 

In  this  case,  S„(t)  =  Sv(t)  =  S(t)  and 

Y  X  ,?0?  (**i 


A    (t)    =      /   S(x)  [cosh(^/ab    (t-x))    -    sinh(-\/ab    (t-x))]dx 
0 

(64) 

t  

Av(t)    =      /  S(t)  [cosh(-\,  ab    (t-x))    -  sinh(  -\/ab"  (t-x) )  ]dx 
Y  0 

(65) 


By  inspection,  A  (t)  >  0  and  A  (t)  >  0  for  t  >  0,  and  both 
biases  are  positive.   This  is  reinforced  by  numerical 
results . 

The  general  conclusion  reached  from  this  analysis  is 
that  the  biases  are  both  positive  unless  one  side  is  going 
to  win  very  decisively.   In  this  case,  the  bias  on  the 
winner's  side  will  be  negative,  while  the  loser's  is  positive 

The  remainder  of  this  section  will  be  a  statement  and 
discussion  of  hypotheses  concerning  bias. 

Hypothesis  2-1 

Given  fixed  initial  force  levels  and  fixed 
attrition  coefficients,  as  the  breakpoint 
force  levels  increase,  the  numerical  bias 
decreases.   But  as  a  percentage  of  casualties 
of  the  deterministic  model,  the  bias  increases. 
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This  hypothesis  has  heuristic  appeal  when  the  bias 
equations  are  analyzed.   The  bias  terms  are  (see  Appendix  D) : 


n 
o 


Sv(t)  =  a[  E   nP(t,m   ,n)  +  n,  m.  P(t,m  +  l,n.  ) 
Y         _         np       Dp  Dp     np    bp 


n=n, 
bp 


mo 

+  n.    E      P(t,m,n,)  ]  (54) 

bP  m=m  +1       bP 
np 


mo 


S  (t)  =  b[  E   mP(t,m,nb  )  +  m  i^  *it,%    ,n   +1) 
m=m,  c  c  c        >. 

n 
o 

+  m.  „   E     P(t,m   ,n)  ]  (55) 

tp  n=nK  +1     np 
bp 

As  n,   and  nv   increase,  the  number  of  terms  in  each  summation 
decreases.   Although  the  second  and  third  terms  are  multi- 
plied by  n.    (which  is  increasing) ,  it  would  seem  that  this 
may  be  more  than  offset  by  the  decrease  in  number  of  terms. 
Even  if  this  is  true,  as  a  percentage,  the  bias  will 

increase  unless  the  actual  bias  decreases  as  fast  as  m  -m,   . 

o   Dp 

Although  parametric  analysis  cannot  logically  prove  the  truth 

of  this  statement,  it  shows  that  the  statement  is  certainly 

plausible . 

As  all  parameters  except  m,   and  n,   are  fixed,  a 

numerical  parametric  analysis  can  be  performed.   A  range 

of  fixed  parameters  was  chosen;  small  and  medium  values  of 

m  and  n  ,  and  several  combinations  of  a  and  b.   In  all 
o      o 

cases  the  results  generally  agreed.   The  following  figures 
illustrate  the  results.   They  plot  A  %  and  A  %  as  defined  in 

A  x 
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Appendix  C  as  functions  of  scaled  time  as  defined  in  the 
same  appendix.   The  important  things  to  notice  are  when  the 
lines  separate  and  how  large  the  final  separation  is  for 
the  various  input  parameters  listed  on  each  figure.   The 
figures  are  summarized  in  Figures  30  and  31. 

In  our  examples,  X  is  taken  always  to  be  the  loser  and 
Y  to  be  the  winner.   As  the  breakpoints  are  increased,  the 
difference  in  X  survivors  decreases  from  slightly  over  three 
to  slightly  over  one.   This  seems  to  confirm  the  first  part 
of  the  hypothesis.   Additionally,  the  rate  of  change  of  the 
bias  is  less  than  the  rate  of  change  of  the  breakpoints,  so 
the  percentage  of  bias  increases,  as  indicated  in  Figure  31. 
This  seems  to  confirm  the  second  part  of  the  hypothesis. 
However,  the  bias  for  the  winning  side  does  not  follow  the 
hypothesis  (Y  wins  very  decisively  in  this  case) .   Indeed, 
it  starts  negative,  becomes  positive,  and  then  decreases. 
But  two  things  must  be  noted  before  the  hypothesis  is  re- 
jected:  First,  the  slope  of  the  AY  curve  is  small;  and 
second,  it  eventually  decreases.   This  leads  to  the  Ay% 
curve  in  Figure  31,  where  AY%  is  increasing. 

It  is  argued  in  Appendix  C  that  A%  is  the  more  important 
measure  of  the  bias.   If  this  is  accepted,  then  the  second 
half  of  the  hypothesis  is  more  important  and,  as  indicated 
above,  it  is  true.   All  of  the  results  in  Table  2  reinforce 
this . 

One  final  way  of  stating  this  conclusion  is;  the  greater 

the  differences  m  -m,   and  n  -n,   ,  the  less  significant  the 

o  bp      o  bp  3 

biases . 
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a 

b 

m 
o 

n 
o 

2b 

Ik 

^£ 

V 

^Y 

v 

0.008 

0.004 

40 

40 

0 

0 

3.22 

8.05 

-0.23 

-  1.96 

0.008 

0.004 

40 

40 

8 

8 

2.62 

8.19 

0.18 

1.61 

0.008 

0.004 

40 

40 

16 

16 

2.11 

8.79 

0.42 

4.40 

0.008 

0.004 

40 

40 

24 

24 

1.64 

10.25 

0.52 

7.41 

0.008 

0.004 

40 

40 

32 

32 

1.15 

14.38 

0.47 

12.43 

0.008 

0.004 

24 

24 

0 

0 

2.48 

10.33 

-0.22 

-  3.13 

0.008 

0.004 

24 

24 

4 

4 

2.10 

10.50 

0.06 

0.88 

0.008 

0.004 

24 

24 

9 

9 

1.68 

11.20 

0.29 

4.93 

0.008 

0.004 

24 

24 

14 

14 

1.31 

13.10 

0.38 

8.74 

0.008 

0.004 

24 

24 

19 

19 

0.94 

18.80 

0.38 

16.10 

0.004 

0.0015 

40 

40 

16 

16 

2.06 

8.58 

0.28 

4.06 

0.004 

0.0015 

40 

40 

24 

24 

1.62 

10.13 

0.36 

7.02 

0.004 

0.0015 

40 

40 

32 

32 

1.12 

14.00 

0.33 

11.79 

0.004 

0.0015 

24 

24 

0 

0 

2.28 

9.50 

0.13 

-  2.58 

0.004 

0.0015 

24 

24 

4 

4 

1.98 

9.90 

0.04 

0.82 

0.004 

0.0015 

24 

24 

9 

9 

1.63 

10.87 

0.19 

4.48 

0.004 

0.0015 

24 

24 

14 

14 

1.28 

12.80 

0.27 

8.49 

0.004 

0.0015 

24 

24 

19 

19 

0.90 

18.00 

0.25 

14.37 

TABLE  II.  Bias  for  Different  Breakpoints 
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Hypothesis  2-2 

For  fixed  F   and  F  ,  and  fixed  attrition 
rate  coefficients,  the  larger  the  initial 
forces,  the  greater  the  numerical  bias, 
but  the  smaller  the  percentage  bias. 


A  generally  accepted  conclusion  for  Lanchester-type 
models  is  a  type  of  "law  of  large  numbers".   For  a  battle 
to  annihilation,  the  differences  between  a  deterministic 
model  and  the  equivalent  stochastic  model  "disappear"  as 
the  initial  force  levels  increase.   Hypothesis  2-2  is  a  more 
general  statement  of  this  not  for  battles  to  annihilation, 
but  to  fixed  casualty/initial  force  level  breakpoints. 
Additionally,  it  is  concerned  only  with  the  differences  in 
force  level  histories.   Hypothesis  1-2  was  a  discussion  of 
the  "law  of  large  numbers"  for  the  probability  of  winning. 

There  is  an  intuitive  appeal  for  the  hypothesis.   The 
time  history  of  the  force  levels  for  a  deterministic  model 
is  something  of  an  "average"  time  history  of  the  stochastic 
model.   As  the  number  of  combatants  involved  increases,  the 
greater  the  aggregation,  and  the  closer  this  average  is  to 
the  mean  force  level  history  of  the  stochastic  model. 

All  but  two  parameters  are  held  constant,  so  a  numerical 
parametric  analysis  can  easily  be  done.   Figures  31  to  35 
show  some  of  the  results  of  the  variation  in  parameters. 
Figures  37  and  38  summarize  them  and  show  that,  for  these 
fixed  parameters,  the  hypothesis  is  true. 

Graphs  of  the  results  in  Table  III  would  give  similar 
results  and  further  reinforce  the  hypothesis. 
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a 

b 

F 
X 

FY 

m 
o 

n 
o 

*X 

V 

Sr 

v 

0.008 

0.004 

1.0 

1.0 

12 

12 

1.75 

14.58 

-0.20 

-  5.70 

0.008 

0.004 

1.0 

1.0 

24 

24 

2.48 

10.33 

-0.22 

-  3.13 

0.008 

0.004 

1.0 

1.0 

40 

40 

3.22 

8.05 

-0.23 

-  1.96 

0.008 

0.004 

0.8 

0.8 

24 

24 

2.10 

10.50 

0.06 

0.88 

0.008 

0.004 

0.8 

0.8 

40 

40 

2.62 

8.19 

0.18 

1.61 

0.008 

0.004 

0.6 

0.6 

24 

24 

1.68 

11.20 

0.29 

4.93 

0.008 

0.004 

0.6 

0.6 

40 

40 

2.11 

8.79 

0.42 

4.40 

0.008 

0.004 

0.4 

0.4 

24 

24 

1.31 

13.10 

0.38 

8.74 

0.008 

0.004 

0.4 

0.4 

40 
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0.52 
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0.008 
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0.2 

0.2 

24 

24 

0.94 

18.80 

0.38 
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1.0 
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0.004 

0.0015 
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24 
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-  2.58 

0.004 
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0.6 

0.6 

24 

24 

1.63 
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0.19 
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0.004 

0.0015 

0.6 

0.6 

40 

40 
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8.58 

0.28 

4.06 

0.004 
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0.4 

0.4 

24 

24 

1.31 

13.10 

0.38 

8.74 

0.004 

0.0015 

0.4 

0.4 

40 

40 

1.64 

10.25 

0.52 

7.41 

0.004 

0.0015 

0.2 

0.2 

24 

24 

0.94 

18.80 

0.38 

16.10 

0.004 

0.0015 

0.2 

0.2 

40 

40 

1.15 

14.38 

0.47 

12.43 

TABLE  III.  Bias  for  Different  Initial  Force  Levels 
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The  general  result  of  this  hypothesis  is  that,  if  the 

forces  are  of  any  appreciable  size  (greater  than  20) ,  with 

breakpoint  force  levels  such  that  m  -m,   and  n  -n,   both 
c  o  bp      o   bp 

exceed  20,  the  deterministic  model  is  probably  an  adequate 
representation  of  the  complex  random  process  of  combat. 


Hypothesis  2-3 

If  any  parameter  is  varied  to  bring  the 
forces  closer  to  parity,  the  time  of  the 
battle  increases  and  the  bias  increases . 


This  hypothesis  is  a  result  of  an  analysis  of  the  solu- 
tion for  Ay(t)  and  Ay(t) .   As  was  shown  previously / 


AY(t)  =   /  [S  (T)cosh(yab(t-T))-S  (x)  "V  a/bsinh  ( -\,  ab(t-T)  ]dT 
0 

(60) 


Ay(t)  =  f    [Sx(x)cosh(  -/Ib(t-T))-SY(T)  -va/bsinh(Vab(t-T)]dT 

(61) 


It  appears  that,  as  t  increases,  the  integral  will  increase 
in  absolute  value.   As  only  one  parameter  is  varied,  numeri- 
cal parametric  analysis  was  easy.   All  results  computed 
indicated  the  hypothesis  is  true.   Figure  40  shows  typical 
results.   It  is  interesting  to  note  that  the  hypothesis 
holds  for  the  loser.   However,  it  is  not  clear  what  happens 
for  the  winner,  but  it  would  appear  that  it  is  a  continuous 
function.   If  it  is,  then  at  some  point,  Ay  reaches  a 
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maximum  negative  value,  and  then  increases  to  a  positive 
value  at  parity.   It  then  appears  the  hypothesis  only  holds 
for  the  loser.   The  results  in  Table  IV  show  the  same 
trends  for  the  winner  and  loser. 

Hypothesis  2-4 

If  a  battle  were  to  be  terminated  at  a 
time  t<0.5tf,  where  tf  is  the  completion 
time  of  the  deterministic  model,  there  is 
no  significant  difference  between  the 
results  of  the  two  models . 

Decisions  to  terminate  battles  are  not  always  made 
because  too  many  casualties  have  been  taken.   A  good  example 
of  this  would  be  the  rear  guard  action  of  a  division  involved 
in  a  retrograde  movement.   The  rear  guard  delays  the  enemy 
force,  in  many  cases,  long  enough  to  allow  the  main  force 
to  retrograde  and  set  up  a  defensive  position,  or  at  least 
allow  a  fresh  rear  guard  to  prepare  for  a  further  delaying 
action.   In  either  case  the  rear  guard  commander  will  be 
tasked  to  delay  the  enemy  for  a  fixed  amount  of  time.   When 
that  time  limit  is  reached,  he  will  disengage,  regardless 
of  whether  he  feels  he  could  engage  the  enemy  longer.   In 
fact,  he  may  have  taken  very  few  casualties.   This  would  be 
a  situation  in  which  hypothesis  2-4  would  be  applicable. 

An  evaluation  of  A  (t)  and  Ay(t)  will  show  that  they 
are  generally  increasing  (in  absolute  terms) ,  with  increasing 
time.   A  study  of  various  combinations  of  parameters  has 
shown  this  to  be  the  case.   The  question  then  is,  at  what 
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bmo2(l-(l-Fx)2) 


a 

b 

0.004 

m 
o 

24 

2b 

19 

n 
o 

24 

2s 

0 

0.88 

0.35 

27 

ano2(l-(l-FY)2) 

0.008 

0.18 

0.008 

0.004 

24 

19 

24 

0 

1.30 

0.37 

58 

0.32 

0.008 

0.004 

24 

9 

24 

0 

1.68 

0.29 

91 

0.42 

0.008 

0.004 

24 

4 

24 

0 

2.10 

0.06 

126 

0.48 

0.008 

0.004 

24 

0 

24 

0 

2.48 

-0.22 

156 

0.50 

0.004 

0.001 

6 

0 

6 

0 

1.06 

-0.07 

275 

0.25 

0.004 

0.001 

12 

0 

6 

0 

3.22 

1.93 

493 

1.00 

0.004 

0.001 

6 

0 

6 

0 

1.06 

-0.07 

275 

0.25 

0.004 

0.0015 

6 

0 

6 

0 

1.13 

-0.10 

291 

0.375 

0.004 

0.002 

6 

0 

6 

0 

1.23 

-0.14 

312 

0.50 

0.004 

0.004 

6 

0 

6 

0 

1.93 

1.93 

oo 

1.00 

0.004 

0.004 

6 

0 

12 

0 

1.06 

-0.11 

137 

0.25 

0.002 

0.004 

6 

0 

12 

0 

1.18 

-0.32 

312 

0.50 

0.0015 

0.004 

6 

0 

12 

0 

1.34 

-0.47 

493 

0.75 

0.001 

0.004 

6 

0 

12 

0 

1.97 

3.22 

00 

1.00 

TABLE  IV.  Bias  for  Different  Attrition  Rate  Coefficients 
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time  into  the  battle  do  A  (t)  and  A  (t)  become  significant? 

A  x 

No  method  was  found  to  determine  this,  but  Figures  41  to  69 

are  typical  of  all  results  found. 

As  might  be  expected,  the  time  at  which  the  two  time 

histories  begin  to  diverge  is  strongly  related  to  the  final 

differences  A  %(tf)  and  Ay%(tf).   As  has  been  noted,  as 

m  -m,   and  n  -n,   decrease,  and  the  force  levels  are  moved 
o  bp      o  bp         ' 

closer  to  parity,  A  %(tf)  and  A  % (tf )  increase.   The  worst 
examples  found  did  not  diverge  before  0.5  tf,  with  one 
excaption.   The  one  exception  to  this  hypothesis  is  when 
the  forces  are  at  or  very  near  parity.   In  this  case,  the 
hypothesis  does  not  hold  because,  with  the  forces  near  or 
at  parity,  the  rate  at  which  casualties  on  both  sides  occur 
approaches  zero  as  time  increases.   This  gives,  in  a  sense, 
an  exaggerated  time  length  of  the  battle  and  thereby  distorts 
the  relative  amount  of  time  Av%  and  A  %  are  appreciable. 

A  X 

C.   VARIANCE 

It  is  possible  for  the  two  models  to  agree  on  which  side 
is  going  to  win,  and  on  the  time  histories  of  the  average 
force  levels,  but  the  deterministic  model  may  not  adequately 
represent  the  complex  random  process  of  combat.   This  is 
because  the  outcome  of  the  stochastic  model  is  determined 
by  the  outcome  of  many  uncertain  events.   The  probability 
distribution  function  related  to  these  uncertainties  may 
cause  the  variance  of  the  force  levels  to  be  considerable. 
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If  the  system  being  evaluated  includes  measures  of  effec- 
tiveness that  are  functions  of  the  distribution  of  the  final 
force  levels,  then  in  the  case  of  large  relative  variance, 
the  deterministic  model  is  not  an  adequate  representation 
of  the  complex  random  process  of  combat. 

Clark  [Ref.  8]  analyzed  the  variance  of  the  square-law 
stochastic  attrition  process  and  hypothesized  there  are 
two  types  of  time  histories  of  the  variance.   The  first  is 
where  the  variance  grows  and  converges  asymptotically  to  a 
final  level,  while  the  second  is  where  it  grows  to  a  peak 
level  and  then  decreases  asymptotically  to  a  final  level 
(see  Figures  70  and  71) . 

The  qualitative  pictures  of  the  time  state  probabili- 
ties in  Chapter  II  provide  some  useful  insights  into  why 
the  variance  asymptotically  approaches  a  value  for  every 
set  of  input  parameters.   As  was  noted,  the  change  in  the 
time  state  probability  is  a  combination  of  diffusion  and 
convective  transport  of  the  probability  mass.   The  diffusion 
provides  the  variance  in  the  force  level.   As  the  probability 
mass  "hits"  the  boundary  it  remains  there  (is  absorbed) ; 
for  the  probability  mass  absorbed,  the  diffusion  effect 
has  also  ended  and  the  variance  of  that  probability  mass 
is  fixed.   As  more  is  absorbed,  the  less  relative  impact  the 
remaining  probability  mass  will  have  on  the  variance.   Even- 
tually, as  time  grows,  the  remaining  probability  will  have 
no  effect,  thus  giving  the  asymptotic  characteristic. 
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However,  because  of  the  convective  transport  of  probability 
mass,  the  remaining  probability  mass  will  not  only  continue 
to  be  diffused,  it  will,  in  general,  arrive  at  the  boundary 
closer  to  the  breakpoint  values  (e.g.,  for  the  X  force,  the 
center  of  mass  of  the  small  probability  mass,  AP, ,  arriving 
in  the  interval  (t,t+At,),  will  be  at  m  survivors,  but 
AP2  arriving  in  the  interval  ( t+At. , t+AtO  will  be  centered 
at  m-k  survivors).   Thus,  as  time  increases,  the  convective 
transport  of  the  probability  mass  causes  the  average  of  the 
probability  mass  not  absorbed  to  approach  the  breakpoint 
force  level.   The  breakpoint  force  level  acts  as  a  boundary, 
and  the  probability  mass  cannot  be  diffused  past  that  point. 
In  fact,  the  probability  that  would  have  been  diffused  past 
the  breakpoint  "piles  up"  at  it,  and  the  variance  is  not 
quite  as  great  as  it  might  be.   In  fact,  if  the  probability 
mass  is  transported  fast  enough  toward  the  breakpoint,  the 
variance  of  the  remaining  probability  mass  may  decrease  in 
spite  of  the  diffusion  and  might  cause  the  overall  variance 
to  decrease. 

For  illustration,  take  Figure  3  and  collapse  the 
probability  mass  to  the  X  (m  survivor)  axis.   The  result 
is : 
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P(t=77.5.m) 


m  survivors 
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Collapsed  Probability  Mass 
FIGURE  72 
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This  gives  a  qualitative  picutre  of  P ( t, 77 . 5 ,m) .   The  vari- 
ance is  calculated  from  the  sum  of  the  squared  deviations 
from  the  mean.   In  this  case,  it  appears  the  variance  will  be 
"large"  (in  fact,  it  is  on  the  order  of  magnitude  of  the  mean) 

If  we  now  take  Figure  5  and  collapse  the  probability 
mass  to  the  X  (m  survivor)  axis,  the  result  is: 


P(t=155.5,m) 


tn  survivors 


iflT 


"3u~ 


C^LLAPSm    pRGnA".U,rrY  MASS 
Figure    73 
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Obviously  the  variance  is  much  small  than  it  was  at 
t  =  77.5. 

Note  that  X  is  the  decisive  loser  in  this  battle;  the 
convective  transport  has  moved  the  probability  mass  toward 
the  breakpoint  fast  enough  to  offset  the  diffusion  effect, 
and  the  variance  has  decreased.   On  the  other  hand,  if  the 
same  graphs  were  made  for  the  Y  force  (the  decisive  winner) 
it  would  be  obvious  that  the  probability  mass  is  not  being 
transported  toward  the  Y  breakpoint  very  rapidly,  and  the 
variance  continues  to  increase. 

Thus  it  seems  that  how  the  variance  progresses  is  depen- 
dent on  how  rapidly  the  side  is  being  attrited,  relative  to 
the  other  side.   If  a  side  wins,  the  variance  will  grow 
asymptotically,  whereas  if  a  side  loses,  the  variance  will 
peak  and  then  decrease  asymptotically.   The  results  of  all 
combinations  of  parameters  in  this  study  indicate  this  type 
of  performance. 

Furthermore,  the  results  in  Table  V  indicated  two 
general  trends . 

Hypothesis  3-1 

As  m  -m,   increases,  the  variance  increases 
in  numerical  value,  but  decreases  in  relative 
measure  to  the  expected  casualties. 


Hypothesis  3-2 

As  the  forces  approach  parity,  the  variance 
increases . 
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Vv(x) 
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0 

4.50 

4.67 

8 

3.66 

4.18 

16 

2.97 

3.61 

24 

2.32 

2.94 

32 

1.71 

2.04 

0 

3.45 

3.59 

4 

2.90 

3.27 

9 

2.34 

2.81 

14 

1.85 

2.26 

19 

1.35 

1.56 

0 

2.39 

2.49 

4 

2.90 

3.27 

9 
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14 

1.82 
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19 

1.33 

1.65 
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2.89 

2.94 

24 

2.29 

2.46 

32 
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1.81 
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3.17 

2.79 
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2.60 

9 

2.28 

2.29 

14 

1.77 

1.92 

19 

1.28 

1.40 

0 

2.20 

1.96 

TABLE   V.      Standard   Deviation   of   Force   Levels 
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The  second  trend  is  in  general  agreement  with  previous 
conclusions  in  the  respect  that,  unless  the  forces  are  near 
parity,  the  models  are  not  significantly  different.   The 
first  trend  is  in  agreement  with  Hypothesis  2-1,  so  that 
the  conclusions  reached  for  it  seem  to  apply  also  for  the 
variance.   Thus,  even  though  the  variance  is  a  different 
measure  for  comparison  of  the  two  models,  it  seems  to  agree 
with  previous  conclusions. 
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IV.   SUMMARY 

The  intent  of  this  analysis  was  to  consider  the  question, 
can  the  complex  random  process  of  combat  be  adequately 
represented  by  a  deterministic  model.   The  models  chosen 
for  comparison  were  the  deterministic  and  stochastic  attri- 
tion versions  of  the  Lanchester  square-law  attrition  pro- 
cess.  Three  aspects  of  these  processes  have  been  compared 
and  analyzed. 

Deterministic  Model  Stochastic  Model 

A.  Fixed  winner  and  loser        Probability  of  winning 

B.  Time  history  of  the  Time  history  of  the 
force  levels  expected  force  levels 

C.  No  variance  of  the  Variance  of  the  force 
force  levels  levels 

Conclusions  in  the  form  of  hypotheses  were  discussed  for 
each  of  the  above  comparisons: 

A.   PROBABILITY  OF  WINNING 

Hypothesis  1-1 

As  the  breakpoint  force  levels,  m,   and 
are  moved  closer  to  the  initial 


force  levels,  m  and  n  ,  respectively, 
the  difference  in  the  probability  of 
winning  between  the  models  increases. 


Hypothesis  1-2 

For  fixed  breakpoint  casualty/initial  force 
ratios   (F   and  Fy) ,  the  difference  between 
the  probability  of  winning  for  the  deter- 
ministic model  and  the  stochastic  model 
decreases  as  the  initial  force  levels  increase 
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Hypothesis  1-3 

As  the  forces  move  away  from  parity,  the 
difference  between  the  probability  of 
winning  for  the  two  models  becomes 
negligible. 


B.   TIME  HISTORY  OF  EXPECTED  FORCE  LEVELS 

Hypothesis  2-1 

Given  fixed  initial  force  levels  and  fixed 
attrition  coefficients,  as  the  breakpoint 
force  levels  increase,  the  numerical  bias 
decreases.   But  as  a  percentage  of  casual- 
ties of  the  deterministic  model,  the  bias 
increases . 


Hypothesis  2-2 

For  fixed  F   and  F„,  and  fixed  attrition 
rate  coefficients,  the  larger  the  initial 
forces,  the  greater  the  numerical  bias, 
but  the  smaller  the  percentage  bias. 


Hypothesis  2-3 

If  any  parameter  is  varied  to  bring  the 
forces  closer  to  parity,  the  time  of  the 
battle  increases  and  the  bias  increases. 


Hypothesis  2-4 

If  a  battle  were  to  be  terminated  at  a  time 
t  _<  0 . 5  tf,  where  tf  is  the  completion  time 
of  the  deterministic  model,  there  is  no 
significant  difference  between  the  results 
of  the  two  models. 

VARIANCE  OF  FORCE  LEVELS 

Hypothesis  3-1 

As  m  -m,   increases,  the  variance  increases 
in  numerical  value,  but  decreases  in  relative 
measure  to  the  expected  casualties. 
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Hypothesis  3-2 

As  the  forces  approach  parity,  the  variance 
increases . 


The  hypotheses  were  formulated  and  evaluated  through 
analyses  of  combinations  of  analytical  solutions,  numerical 
solutions,  and  graphical  representations. 
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V.   DISCUSSION 

The  hypotheses  formulated  in  this  thesis  provide  a  basis 
for  the  extension  of  several  conclusions  about  the  adequacy 
of  the  deterministic  model.   When  one  side  is  going  to  win 
decisively  in  the  deterministic  model,  the  models  do  not 
differ  significantly  and  the  deterministic  model  is  to  be 
preferred.   How  decisively  one  side  must  win  to  make  the 
difference  insignificant  is  a  function  of  two  sets  of 
parameters;  the  initial  force  levels  and  the  breakpoint 
force  levels.   As  the  initial  force  levels  increase,  the 
relative  differences  between  models  decrease  as  long  as 
m  -m,   and  n  ~n,   do  not  increase.   Thus,  with  relatively 
large  initial  force  levels  (greater  than  20)  and  breakpoint 
force  levels  such  that  there  will  be  relatively  large  losses 
allowed  (also  greater  than  20) ,  the  models  are  not  signifi- 
cantly different. 

It  must  be  recognized  that  these  conclusions  only  apply 
to  the  very  simple  and  idealistic  Lanchester  square-law 
attrition  process  and  the  square*law  stochastic  attrition 
process.   As  is  always  the  case,  the  insights  into  combat 
dynamics  obtained  from  these  models  are  no  more  valid  than 
the  models  themselves.   However,  intuition  indicates  the 
conclusions  may  be  generalized.   It  is  intuitively  appealing 
to  say  that,  no  matter  what  the  model  or  degree  of  complexity, 
if  the  deterministic  model  predicts  that  one  side  will  win 
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decisively,  the  stochastic  model  will  not  disagree 
significantly . 

Nevertheless,  it  is  felt  that  the  hypotheses  may 
provide  a  point  of  departure  for  the  comparison  of  more 
complex  models  of  combat. 
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VI.   CONCLUSIONS 

For  the  idealized  model  (i.e.,  the  so-called  Lanchester 
"square-law"  attrition  process) ,  it  is  concluded  that  the 
complex  random  process  of  combat  can  be  adequately  repre- 
sented by  a  deterministic  model  if  the  following  conditions 
are  met: 

1.  Each  side  starts  with  at  least  20  combatants. 

2.  Each  side  is  willing  to  take  at  least  20  casualties. 

3.  The  forces  are  not  near  parity. 

4.  If  the  forces  are  near  parity,  but  each  side 
initially  has  in  excess  of  40  combatants  and 
is  willing  to  take  in  excess  of  20  casualties. 

Consequently,  it  seems  plausible  that  similar  results  would 
be  obtained  if  one  were  to  compare  analogous  deterministic 
and  stochastic  "real  world"  combat  models  (i.e.,  large- 
scale  complex  system  models  of  military  forces) .   In  other 
words,  it  seems  reasonable  to  expect  that  if  the  opposing 
military  forces  are  not  near  parity  and  have  relatively 
"large"  numbers  of  initial  combatants,  a  deterministic  model 
may  be  used  to  model  the  dynamics  of  combat  without  losing 
any  essential  information  about  the  battle's  outcome.   From 
the  research  reported  here,  the  author  concludes  that  one 
may  safely  take  "large"  numbers  of  initial  forces  (in  par- 
ticular, of  tanks) ,  to  mean  units  of  battalion  size  or 
larger.   Thus,  it  appears  as  though  current  U.S.  Army  deter- 
ministic models  such  as  the  Bonder/IUA  model  are  to  be 
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preferred  for  sensitivity  analysis  because  of  their  compu- 
tational advantages,  once  they  are  "calibrated"  from  Monte 
Carlo  simulation  output  (for  example,  from  DYNTACS  or 
CARMONETTE) . 
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APPENDIX  A 
COMPARISON  OF  TWO  METHODS  OF  OBTAINING  THE  PROBABILITY  OF  WINNING 

In  this  appendix  it  is  shown  that  two  apparently  differ- 
ent methods  of  obtaining  the  win  probabilities  are  really 
the  same.   In  his  doctoral  thesis,  Springall  [Ref.  28] 
derived  a  recursive  expression  for  the  probability  of 
winning  in  a  homogeneous  force,  constant  coefficient  sto- 
chastic Lanchester  model.   His  expressions  and  the  develop- 
ment of  them  are  quite  different  from  those  of  Brown  [Ref.  7]. 
Brown  developed  the  familiar  set  of  partial  difference 
equations : 


P(m,n:r,s)  =  A(r,  s'^B  (r,s)  P  (■"."'*-!, s) 


+  A(r,s)'r;SB(r,s)  P(m,~:r,B-l)    (66) 


with  boundary  conditions 


1    n  =  s 

P(m.  ,n:m,  ,s)  =  (  (67) 

bp    bp      ^0    otherwise 


1   m  =  r 
10    otherwise 


P(m'nbP,r'nbp>  =  L    _^._J„  (68) 


P(nv  ,n:r,s)  =0   n  >  s;     P(m,n,   :r,s)  =0   m  >  r 

(69) 
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1   m,   =  r   and  n,  _  =  s 

vnbP:r's)  \n   :i  .  (7o) 

otherwise 


where  P(m,n,:m  ,n  )  =  probability  X  will  win  with  m 
survivors.   When  this  is  summed  from  m.  +1  to  m  ,  it  yields 
the  probability  X  will  win: 

m 

Px  =   Z      P(m,n   :mo,no)  (71) 

m=m,  +l 

When  the  above  equations  are  looked  at  carefully,  what 
they  are  doing  is  explained  by  the  following  figure  (Fig. 
74)  . 

The  probability  that  X  will  win  with  m  survivors  is 
simply  the  sum  of  the  probabilities  of  all  possible  paths 

from  (m  ,n  )  to  (m,n,  )  with  the  probability  of  going  from 

A  ( r  s) 
state  (r,s)  to  state  (r-l,s)  equal  to  A(r  s\  +  s(r  s)  and 

the  probability  of  going  from  state  (r,s)  to  state  (r,s-l) 

equal  to  -r-7 1     ' SL  , x  .   And  the  probability  X  will  win 

A^r,sj  +  b(r,s; 

is  simply  the  probability  that  X  will  win  with  m  survivors, 
or  m  -1  survivors,  or  ...,  or  m,  +1  survivors. 

A  similar  analysis  of  Springall's  equations  will  show 
that  the  two  methods  are  identical.   Springall  concludes 
tli  at 


P(m,nfa  :mQ,no)  =  B  (m,!^  +1)  X  (m,!^  +1)  (72) 
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Method  of  Solving  Recursive   Relationship  for  the   Probability 
of  Winning   for  Brown's   Equations 

FIGURE   74 
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where  X(m,n,  +1)  is  derived  from  the  following  set  of 
recursive  relations; 


X(mo'no}  =  A(m  ,n  )  +  B(m  ,n  )  (73) 

o  o       o  o 


B(m  /n+l)X(m  ,n+l) 

X(m  ,n)   =  _  .  ° — ,  ,  _  ,° r —    n,   <  n  <  n      (74) 

o       A(m  ,n)  +  B(m  ,n)       bp        o      v 


A(m+l,n  )X(m+l,n  ) 

X(m,n  )   =  •=-! r — ; — =-7 r —    m    <  m  <  m       (75) 

'  o     A(m,n  )  +  B(m,n)      np        o 


Y,    y  A(m+l,n)X(m+l,n)  +  B (m,n+l) X (m,n+l)     ,-,. 

Aviu.n;    —  r—3 r — : — =t"7 ^ WO/ 

'  A(m,n)  +  B(m,n) 


Define 


P'(r,s)  =  (A(r,s)  +  B  (r  ,s)  )  X(r  ,s)  (77) 


Then 


P,(mo,n0)  =  1  (78) 


B(m  ,n+l) 

P'  (m  ,n)   =  ^-t ~n  t  5-7 -rr  P1  (m  ,n+l)       (79) 

o       A(m  ,n+l)  +  B(m  ,n+l)      o' 


A(m+l,n  ) 
p,^V   s  A(m+l,n  )  +B(m+l,n  l^'^V       (80) 
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n,, x  _  A(m+l,n)P'  (m+l,n)    ,  B  (m,  n+1)  P  '  (m,  n+1) 


A(m+l,n)  +  B(m+l,n)    A(m,n+1)  +  B(m,n+1) 

(81) 


Therefore 


A(m.  +  l,n)P'  (iru  +l,n) 

P(m,  ,n:m  ,n  )  =  -  .    bP  ■ >  t  d  ,P — —, r  (82) 

top'    o'  o    A (n^  +1,11)  +  B(itl  +l,n) 


A(m,n,  +  l)P*(m,n,  +1) 

P(m,n.  n:m,n  )  =      bp,       ,  bP  (83) 

bp  o  o    A(m,n   +1)  +  B(m,n.  +1) 


What  these  equations  say  is,  the  probability  of  starting 

with  (m  ,  n  )  survivors  is  one.   The  probability  of  having 

A  ( m  ,  n  ) 
o   o 

one  more  casualty  from  X  is  ■=—> c — ; — sn \   and  from 

tw  x  A^m    ,n ) +   B(m    ,n    ) 

B(m    ,n    )  o      o  o'    o 

Y    is      —, r — ; — ^-7 .  .      The   probability   of   having    (m,n) 

A(m  ,n  )  +  B(m  ,  n  )        c  2  r 

o  o       o  o 

survivors  is  the  probability  of  having  (m+l,n)  survivors 
times  the  probability  of  an  additional  casualty  being  from 
X,  added  to  the  probability  of  having  (m,n+l)  survivors 
times  the  probability  of  an  additional  casualty  being  from 
Y.   This  recursive  relationship  leads  to  the  following 
figure  (Figure  75) ,  which  is  the  same  as  Figure  74.   Thus, 
even  though  the  sets  of  relationships  appear  to  differ, 
they  do  yield  the  same  results.   This  was  further  supported 
by  a  comparison  of  results  from  Lee  and  Wannasilpa  [Ref .  20] 
and  the  results  from  the  enclosed  computer  program,  which 
was  set  up  using  the  recursive  equations  of  Springall.   For 
the  same  parameters,  the  results  were  the  same  within  0.00001 
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FIGURE  75 

The  differences  are  assumed  to  have  occurred  because  of 
roundoff  errors  within  the  computer  system  used. 
The  program  accuracy  was  further  checked  by 


m 


m=mb: 


n 


P(m,n.   :m  ,n_)  +  I             P(m,   ,n:m  ,n  ) 

,  T       Dp   O   O  x  ftp'     o   O 

+  1     r  n=n,  +  1    ^    ww 

>  bp 


=  1 


(84) 


All  results  were  in  the  interval  (0.9999692,1.0000000). 
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APPENDIX  B 
NUMERICAL  SOLUTION  OF  THE  FORWARD  KOLMOGOROV  EQUATIONS 

Exact  solutions  for  the  time  state  probabilities  for 
the  square-law  stochastic  attrition  process  can  be  developed 
only  for  specific  initial  force  levels.   To  the  author's  best 
knowledge,  no  solution  for  the  time  state  probabilities  for 
general  initial  force  levels  has  been  developed.   Even  if 
one  were  able  to  develop  such  a  general  solution,  the  result 
would  be  of  no  practical  use  due  to  its  complexity.   The 
solution  for  the  linear-law  stochastic  attrition  process 
was  developed  by  Clark  [Ref.  8],  but  it  is  too  complex  for 
practical  use:   So,  solutions,  if  they  could  be  developed, 
for  more  complex  Lanchester-type  stochastic  attrition  pro- 
cesses, would  most  likely  be  of  no  practical  use.   Therefore, 
the  most  practical  way  to  solve  such  equations  is  through 
numerical  integration  techniques.   The  method  chosen  for 
this  thesis  was  the  fourth  order  Runga-Kutta  method  [Ref.  21] , 
The  computer  program  is  set  up  to  solve  Lanchester-type 
attrition  process  models  of  the  form 


^dt™'1^  =  A(m+1'n)p(t'm+1'n)  +  B(m,n+l)P(t,m,n+l) 


-  (A(m,n)  +  B(m,n) )P(t,m,n)        (85) 


14  8 


with 


P(0,m,n)   =  < 


1   n  =  n    and   m  =  m 
o  o 

0    otherwise 


(86) 


P(t,m,n) 


0   m<m    or   n<n 
o  o 


(87) 


and  boundary  conditions 


dP(t,m,rv    ) 

gt  =   B(m/nbp+l)P(t/m,nbp+l) 


(88) 


dP(t,m      ,n) 

ZT2 =  A(mbpM/n)P(t/mb  +  l,n) 


(89) 


However,  with  a  few  minor  alterations,  the  program  will 
solve  the  more  general  Lanchester-type  stochastic  attrition 
processes  with  attrition  rates  of  the  form  F(t,m,n)  and 
G(t,m,n) . 

Numerical  precision  of  the  Runga-Kutta  methods  is 
dependent  on  the  length  of  the  time  step  chosen.   For  the 
time  steps  chosen,  numerical  precision  appeared  to  be 
satisfactory  since 

m     n 
o     o 

E     E    P(t,m,n)  was  always  in  the  interval 

m=nv   n=n, 
op    bp 

(0.9992,1.0003) . 
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Because  the  expected  number  of  survivors  was  calculated 
from 


m  n 
o     o 

m(t)  =  I  I    mP(t,m,n)  (90) 

m=m,  n=n, 
tip    bp 


it  was  felt  the  solutions  were  accurate  enough.   Additionally, 
because  exact  solutions  exist  for  general  initial  force 
levels  for  P(t,m  ,n)  and  P(t,m,n  ) ,  they  were  used  to 
check  the  accuracy  of  the  approximate  solutions .   The 
approximate  solutions  never  deviated  from  the  exact  solution 
by  more  than  0.00005. 
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APPENDIX  C 


SCALING 


In  a  comparison  of  a  stochastic  model  and  the  equiva- 
lent deterministic  model,  two  classes  of  comparisons  need 
to  be  considered.   The  first  is  the  difference  in  results 
obtained  from  the  two  models.   The  second  is  to  take  these 
differences  and  see  how  they  change  with  variations  in 
parameters.   One  runs  into  difficulties  when  trying  to  make 
sense  out  of  these  comparisons  when  results  are  expressed 
in  dimensional  units  (i.e.,  What  is  the  significance  of 
the  fact  that  at  150  minutes  into  a  battle,  the  difference 
in  survivors  is  1.8  tanks?). 

One  method  of  overcoming  this  problem  is  the  technique 
of  scaling.   If  comparisons  are  made  when  results  are  in 
dimensionless  quantities,  it  is  much  easier  to  see  if  there 
are  any  significant  differences.   For  example,  one  can 
easily  compare  two  probability  statements  and  note  any 
difference  between  them. 

Scaling  is  a  subset  of  a  "science"  called  dimensional 
analysis,  a  technique  often  used  in  the  physical  sciences. 
For  example,  in  the  area  of  fluid  dynamics,  Reynolds  deduced 
from  some  experiments  that  laminar  flow  broke  down  into 
turbulent  flow  at  some  critical  velocity  above  that  at 
which  turbulent  flow  was  restored  to  the  laminar  condition; 
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the  former  velocity  being  called  the  upper  critical  velocity/ 
and  the  latter,  the  lower  critical  velocity.   Reynolds  was 
able  to  .generalize  his  conclusions  by  the  introduction  of 
a  dimension less  term  called  the  Reynold's  number.   He  found 
that  certain  critical  values  of  the  Reynold's  number  defined 
the  upper  and  lower  critical  velocities  for  all  fluids 
flowing  in  all  sizes  of  pipes,  and  thus  deduced  that  single 
numbers  define  the  limits  of  laminar  and  turbulent  pipe 
flow  for  all  fluids.   Many  other  dimensionless  numbers  are 
used  to  generalize  results  —  the  Froude  number,  the  Cauchy 
number,  the  Mach  number,  etc.   Even  in  mathematics  dimen- 
sionless quantities  are  used  to  generalize  results  —  it  and 
e,  for  example.   (it  is  the  ratio  between  the  circumference 
of  a  circle  and  its  diameter,  regardless  of  the  units  of 
measure  employed.) 

Scaling  involves  the  ratio  of  two  quantities  with  the 
same  dimensions.   The  dimension  of  a  ratio  of  two  quantities 
with  the  same  dimensions  is  said  to  be  dimensionless.   The 
analysis  in  this  paper  is  concerned  with  two  classes  of 
dimensions;  time  and  quantity  (number  of  casualties/survivors, 
and  the  differences  in  casualties/survivors) .   Either  model 
(the  stochastic  or  deterministic)  could  have  been  chosen  as 
the  basis  for  forming  ratios  for  the  two  classes  of  dimensions; 
the  deterministic  model  was  used  as  a  basis.   Time  was  trans- 
formed into  a  dimensionless  quantity  by  using  the  ratio  t/tf, 
where  tf  is  the  time  at  which  the  deterministic  battle  would 
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end.   It  is  easy  to  see  that  if  0  <  t  <  tf ,  then  the  ratio 
is  between  zero  and  one. 

The  transformation  of  the  quantities  requires  some 
explanation  of  the  rationale  behind  it.   The  transformation 

used  was 

.         ,    number  of  casualties  at  t 

ratio       number  of  casualties  at  tf  for  the  deterministic  mode 

It  is  easily  seen  that,  for  the  deterministic  model  with 

0  <  t  <  t^;   0  <  A   . .   <  1.   A  dimensionless  quantity  was 
—   —  r     —  ratio  —  -a      j 

desired  that  would  express  the  "importance"  of  the  number 

of  casualties.   A  difference  of  one  casualty  out  of  forty 

is  not  as  "important"  as  a  difference  of  one  casualty  out 

of  four.   A    .   expresses  this  "importance".   For  example, 

at  tr,  A   ,  .   =0  for  the  deterministic  model,  while 
f   ratio  ' 

A   .  .   ^  0 ,  in  general,  for  the  stochastic  model.   In  fact, 


A  the  difference  in  casualties  at  t£ 

A    .   =   f 

number  of  casualties  at  tf  for  the  deterministic  model 


as  the  number  of  casualties  at  tf  for  the  deterministic 

model  increases,  A   . .   decreases,  giving  the  desired  result. 

'   ratio  '     3 
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APPENDIX  D 

DERIVATION  OF  EXPRESSION  FOR  BIAS  IN  THE  SQUARE  LAW  STOCHASTIC 
ATTRITION  PROCESS  WITH  FIXED  BREAKPOINT  FORCE  LEVELS 


From  previous  results  it  is  known  that  the  square-law 

stochastic  attrition  process  is  defined  by: 

for   m,  +1  <  m  <  m     and  n,  +1  <  n  <  n  ; 

r>p   —   —  o  bp   —   —  o 


dP(tdt'n)  =  anp(t'm+1'n)  +  bmP(t,m,n+l)  -  (an+bm) P ( t,m,n) 

(91) 

and  on  the  boundaries  of  the  state  space, 


dP  ( t ,  m ,  n,  ) 

dt  bp  =  bmP(t,m,nbp+l)  (92) 


dP(t,m   ,n) 

^2 =  anPU^p+^n)  (93) 


also, 


P(t,m,n)  =  0     for  m  >  m    and/or   n  >  n 
v  '  '  o  o 


and 


1    m  =  m   and   n  =  n 
P(0,m,n)  = 

0    otherwise 
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Define 


A  P(t,m,n)  =  P(t,m+l,n)  -  P(t/m/n) 


and 


A  P(t,m,n)  =  P(t,m,n+1)  -  P(t/m/n) 


so  that  (91)  may  be  written  as 


for  m,  +1  <  m  <  m    and   n,  +1  <  n  <  n 
bp   —   —  o  bp   —   —  o 


dP(dtm/n)  =  anAmp(t'm'n)  =  bmAnP(t,mfn)  (94) 


First  consider  two  jointly  distributed  discrete  random 
variables,  M  and  N  with  probability  mass  denoted  as 


f  .-(ni/n)  .   Further,  assume  that  f  N(m,n)  =  0  for  m  <  m, 

and  n  <  n,  .   Th< 
bp 

E(M)  is  given  by 


and  n  <  n,  .   Then  the  expected  value  of  M,  denoted  as 


00  00 


E(M)  =   I      E     mfM  M(m,n) 

m=m,   n=n, 
op    bp 


This  may  be  expressed  in  terms  of  the  conditional  expectation 
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Z     mfM  M(m'n) 
m=m, 

E(M:N=n)  =  2E. 


fN(n) 


as  E(M)  =   Z     f  (n)E(M:N=n) 

n=n, 
bp 


Note  that 


m 
o    n 

§£  =  [ft  E     Z   mP(t,m,n)]  (95) 

m=m,   n=n. 
Dp    bp 


To  compute  -rr,    which  is  given  by  equation  (101)  below,  first 
sum  equation  (9  4)  over  all  n,  from  n,  +1  to  n   to  obtain 

n  n 

d     °  ° 

■g£    [    Z     P(t,m,n)]  =  a   Z      nAmP(t,m,n) 

n=n,  +1  n=n.  +1 

bp  bp 

n 
o 

+  bm  Z      A  P(t,m,n)      (96) 

n=n,  +1 
bp 

Note  that 


o  o 

Z       A  P(t,m,n)  =   Z       [P(t,m,n+1)  -  P(t,m,n)] 

n=n,  +1  n=n,  +1 

bp  bp 


=  P(t,mfno+l)  -  P(t,m,n   +1) 


or 
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n 
o 

E  A   P(t,m,n)    =   -   P(t,m,n,     +1) 

n=n^    +1  n  Dp 


bp 
Thus 


n  n 

-         o  o 

%rr[   E  P(t,m,n)]  =   a[    E  nA  P(t,m,n)] 

UL  -  .    ■  III 

n=n,     +1  n=n,     +1 

bp  bp 


-  bmP(t,m,nbp+l)  (97) 


Adding    (92)    to    (97)    results    in 

n 

j  o 

|^[    E  P(t,m,n)]    =   a[    E  nA   P(t,m,n)]  (98) 

n=n.  n=n,     +1 

Dp  bp 

Multiply  (9  8)  by  m  and  sum  over  m  from  ni  +1  to  m  ,  and 
add  (9  3)  multiplied  by  m,   : 

m     n 

'■Li  E      E     mP(t,m,n)]  =  — 

m=m,   n=n, 
op    bp 

m       n  n 

o       o  o 

=  a[  E        E  mnA  P(t,m,n)]  +  m.    E    anP(t,m,  +  l,n) 
,  -,                 .-,  m  op  op 

m=m,  +1  n=n,  +1  r   n=n,  c 

op      bp  bp 

(99) 


Observe  that 
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m 


n 


m=mbp+1     n=nbp+1 


mnA    P(t,m,n) 


m 


m 


=      E 


n 


mA^      E 
m 


nP  (t ,m,n) 


m=m,     +1  l  n=n,     +1 

op  bp 


m  n  n 

o  o  o 

Z  m[    E  nP(t,m+l,n)    -      E  nP(t,m,n)] 

n=n. +1 

bp 


m=m,     +1      n=n,     +1 


m  n 

o  o 


=      E 


m  E  nP(t,m+l,n) 

n 
bp 


m 


-      E 


n 


m=m,     +1   n=n,__. 


m   E  nP(t,m,n) 

n 
bp 


m=m,     +1  n=n. 


m  +1 
o 


n 


E  (m'-l)       E  nP(t,m',n) 

m=nbp+2 


n=n,     +1 
bp 


m 


n 


o  o 

E  m      E  nP(t,m,n) 

m=nv    +1      n=n,     +1 
b»p  bp 


m  +1 
o 


n 


E  (m'-l)       E  nP(t,m*,n) 

m=rnbP+1  n=nbP+1 


m  n 

o  o 


E  m      E  nP  (11,111,11) 

m=m,     +1      n=n,     +1 
bp  bp 


n 


m,  E  nP(t,m     +l,n) 

*p  n=nK   +1  np 

bp 
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m  n 

o  o 


=   E      (m'-l)   E     nP(t,m',n) 

m=m,  +1       n=n,  +1 
op  op 


m  n 

o  o 

E  m   E     nP(t,m,n) 

m=m,  +1   n=n,  +1 


:mbP+1   n=nbp 


n 
o 


nP(t  ,111  +  l,n) 


v   n=n,  +1      "* 
bp 


m  n 

o  o 


=   z      (m'-l-m)  l  nP(t,m,n) 

m=mbp+1         n=nbp+1 


n 
o 


-  m    z    nP(t,m  +  l,n) 
np    +1      np 

bp 


o      "o 

=   E        E      nP(t,m,n) 

m=m,  +1  n=n,  +1 
bp      bp 

n 
o 


ni    E     nP(t,ni  +l,n)         (100) 


n=n,  +1 
bp 


Substituting  (100)  into  (99) 


m       n 

dm  =      °       ° 
dt 


E 

E             nP(t,m,n) 

1=mbp+1 

n=nbP+1 

n 
o 

-  ^p 

E 
n=nbp+1 

no 

nPft,!^  +l,n) 


+  am.    E    nP(t,m  +  l,n) 

tp  n=nK       nP 
bp 
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nP(t,m,n)    +   arn^  P(t,itL    +l,n     ) 

+  1    n=n,      +1  tr       tr  f  tr 


m 
o 

n 
o 

a      Z 

Z 

m=m,     +1 

n=n, 
bp 

m 
o 

n 
o 

m 
o 

■a      Z  Z        nP(t,m,n)    +   a  Z  n,     P(t,m,n      ) 

m=mbp+l  n=nbp  m=mbp+1     ?  P 


+   arabpnbpP(t'mbP+1'nbP) 


m  n  n 

o  o  ° 

■a      Z  Z        nP(t,m,n)  +   a  Z  nP(t,m,     ,n) 

m=mbp  n=nbp  n=nbp  P 


m 
o 

+  a  E        J.1nbpP(tfm/nbp) 
m=rnbP+1 


+   anbpnbPP(t/mbP+1/nbP) 


n 


dm 


gj.  =   -an  +   a[    Z        nP(t,mb    ,n)    +  n b  m     Plt,^  +l,n      ) 


n=n 


bp 


m 


+   V      Z  P<t,m,nbp>] 

m=rT1bp+1 


(101) 
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Equivalently , 


m 


>P 


m=rnbi 


n. 


r  n=n,     +1  * 

bp 


(102) 
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c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 


COMPUTER    PROGRAM 


SOLUTION    TO    STOCHASTIC    LANCHESTER    EQUATIONS 
USING    FOURTH    ORDER    RUNGA-KUTTA    APPROXIMATIONS 


LO=NUM8 

SOLU 

MC=NUMB 

NC=NUMB 

H=SIZE 

A(X,Y)= 

B(X,Y)= 

FX=BREA 

FY=BREA 

MdP=(l- 

NBF=( 1- 

DI  MENS  I 

DIMENSI 

OIMENSI 

DIMENSI 

DIMENSI 

DIMENSI 

IF    A    PL 

SET 

WHIC 

PR06 

(NOT 


ER    OF    T 
TICN    TD 
ER    OF    C 
ER    OF    C 
OF    TIME 
• «ATTRI 
••ATTRI 
K POINT 
KPOINT 
FX.)*MG 
FYJ*NO 
ON    KX(N 
PROB 
WK(M 
D(L0 
P1IL 
EST( 
OF    T 
IT  =  0. 
H    A    THR 
ABILITY 
£-    TIME 


IME    INTERVALS    NEEDED    FOR    COMPLETE 
DETERMINISTIC    EQUIVALENT     EQUATIO.M 

OMBATANTS    X     STARTS    WITH 

OMBATANTS    Y     STARTS    WITH 
INCREMENT 

TIQN"    RATE    FOR    X     IN    TERMS    OF    X-l    6 

TION"1     RATE    FOR    Y     IN    TERMS    OF    X-l    & 

FORCE    RATIO    FOR    X 

FORCE    RATIO    FOR    Y 


Y-l 
Y-l 


ON 
ON 
GN 
ON 
ON 
OT 


DIMENSICN    PI (313t41),P2(313,41)»EST(313) »D(313) , 

IF  (Zi  ,SIZ£(i2)  ,KX(  100),  KYI  100  J, 

2PPCB(41r4i)tWK(41,41,3), 

3VAR(3i3J ,DETERM(313) 
INTEGER    X,Y»U,Z 

REAL*8    TTL(12)/»CRAIG    J.f,5*'  •  »  ■  SMC    1462,,5*« 

REAL    XTITLE/'TIME'/fYTITLE/'SURV'/iZTITLE/4  «/ 

REAL    KltK2fK3tK4»KlA,K2A,K3AfK4AfK2BtK3BtK4B,KlCfK2C, 

1K1B,K10,K2G,K;>D,K4D,K1E,K2E,K;>E,K4E,K1G,K2G,K3G,K4G, 

2K4F,K3C,K4C,K1F,K2F,K3F 


41) 


C 

C 

LGGICAL*1  IDN(MO 

C 

LCGICAL*1  IDN(41 

F(1)=0.0 

F(2)=0.0 

SIZE( 1)=7.0 

SIZE(2)=9*0 

LC=312 

MC  =  40 

N0  =  40 

NKXY=100 

MC1=M0+1 

LC1=L0+1 

NC1=N0+1 

IT=10 

Al=0.008 

61=0.004 

H=2.0 

FX=1.0 

FY=1.0 

MBF=( 1-FX)*MG 

NBP=(1-FYJ*N0 

MS=M0+1-MBP 

NS=N0+1-NBP 

M6P1=M3P<-1 

NBP1=N8P+1 

MS1=MS-1 

600  DO  13  L=1,L01 
EST(L)=0.0 

162 


VAR(LI=0.0 
13  CONTINUE 
DC  400  N=l 
Y=NO+2-N 
DO  ^OQ  M=l 
X=M0+2-M 
DO  210  L=l 
IML.EQ.1J 
IF(X.EQ.MO 
IF (X. Eg. MB 
1F(Y.EQ.N0 
IFiy.EQ.N3 
K1=A(X*1,Y 

H-B(XtY) )*P 
K2=A(X+1,Y 

1CP2CL|XI«-P 

2*0.5) 
K3=A(X+ 

UP2(L,X 

2*0. 5) 
K4=A(X+ 

1*<P1(L- 
PKLtX) 
GC    TO    2 

500  IF(X.EQ 
P1(L,X) 
GC    TO    2 

501  F1(L,X) 
GO    TO    2 

502  IFiY.NE 
P1(L,X) 
GC    TO    2 

503  IF(Y.EQ 
K1A=B(X 
K2A=B(X 

1-(A(X,Y 
K3A=B(X 

1-(A(X,Y 
K4A=B(X 

1-(A(X,Y 
P1(L,X) 
IF(Y.EG 
GC    TO    2 

601  E=ABS(P 
1 (EXP< (A 

IF(E.LE 
H=H*0.5 
GO    TO    6 

504  IF(Y.EQ 
K1F 

K2F=A(X 
K3F=A(X 
K4F=A(X 
P1(L,X) 
GC    TO    2 

505  PKLtX) 
GC    TO    2 

506  K1D=A( X 
K2D=A( X 

1 (A(X,Y) 
K3D=A( X 

1-(A(X,Y 
K4D=Al X 

1-(A(X,Y 
PKLtX) 
IFiX.EQ 
GO    TO    2 

602  E=ABS(P 
1 (EXPU 6 

IF(E.LE 
H  =  H*-0.5 
GO    TO    6 


,NS 
fMS 

»L01 

GO    TO 
+  1  )G0 
P+DGO 
♦  DGO 
P  +  DGO 
HPilL 
KL-1, 
)*0.5* 
2(L-1, 


1,Y)*0.5* 
)+P2(L-li 

l,Y)*Pl(L 

1,XDH*K3 

=P1(L-1,X 

00 

.MO+l.AND 

=  0.0 

00 

=  1*0 

00 

.NO+DGO 

=  EXP(-(A< 

00 

.NBP+DGO 

,Y*1)*P2( 

,YM)*0.5 

)+B(X,YD 

,YJ-1)*0.5 

J+acXfYi ) 

iY#U*P2( 
)+B<X,YD 

=-P  1  (  L- 1  ,  X 

.NO.ANO.L 
00 

1(2,M01)- 
(MO+ltNO* 

.COODGO 


500 
TO    502 

TO    504 
TO    506 

TO    507 
-1,X  +  D+B(X,Y+D*P2(L-1,X)-<A(X,Y) 
X) 

(Pl(LtX+l)+Pl(L-lflX+l)}+B(XfY+l)*0.5* 
X) |-(A(X,Yf+B(XfY)J*(Pl(L-lfXJ*H*Kl 

(PKLiXfl)  «-PKL-l,X  +  l)  DB(X,Y  +  1)*0.5* 
X) i-(A(Xf Y)*b(X,Y))*(Pi(L-l,X)+H*K2 

iX*-ll+B(XtY*l)*P2(L,X)-<A(X»Y)+B(XtY)J 

) 

)*(H/6.0)*(KD-2.0*K2  +  2.0*K3  +  K4) 

.Y.EQ. NO+DGO    TO    501 


TO    503 
X,Y)+B(X,Y))*(L-1)*H) 

TO    509 
L-lt  X*-(A(X,Y)«-B(XtY))*Pl(L-ltX) 
*{P2(LfX)+P2(L-ltX)J 

*( P1(L-1,X)  *H*K1A*0.5) 
*(P2(L,XJ*P2(L-ltXJJ 

*( P1(L-1»XDH*K2A*0.5J 

L,X) 

*<  P1(L-1,X) +H*K3A) 

)  +  (H/t>.0J*(KlA  +  2.0*K2A+2.0*K3A  +  K4A) 

.EQ.2)G0    TO    601 

(  (  B(MO  +  lfNO*-l)*NO/A<MOt-l,NO  +  l)  )* 
1)/N0)*H)-1.0)*P2( 2»M01) ) ) 
TO    200 


00 
.NBP+ 

+  1,Y) 
+  D  Y) 
+  1,Y) 
=  PHL 

00 
=  0.. 
00 

+lf 

+  1. 

+  BI 
+  D  Y) 
14-BU 
D,Y) 
)+B(X 
=  P1(L 
.MO.  A 
00 
1<2-,M 


DGO    TO    505 

=  AiX*l,Y)*PUL-l,X+l) 
*0.5*(PUL,XD)  +  P1(L-1,X  +  U  ) 
*0.5*(PKLfX*-ll+Pl(L-ltX+i)) 
-Pit  Li X+D 
-ltX)+(H/6.0)*(KlF+2.G*K2F+2.0*K3F+K4F) 


0 

Y) 

Y) 
X, 


(MO  +  1 
.0*00 

00 


*Pl(L-l»X«-lJ-(A(X,Y)+B<X,Y))*Pl(L-ltX) 

*0.5*(Pl(L,X+l)+Pl(L-lfXt-l))- 

Y)  )*<PHL-lfXi+H*KlD*0.5] 

*0.5*(  Pl(LiX«-n+Pi(L-l  ,X*D  ) 

,Y))*(PHL-1,X) +H*K2O*0.5) 

*P1(L»X+ 1) 

,Y)) *( PHL-lfX)+H*K3D) 

-liX  )  +  (H/6.0)*(KlD+2.Q*K2D  +  2.0*K3D+K4D) 

ND.L .EC.2)  GO    TO    602 

0)-( (  A(MO  +  1,NO  +  D*MO/B(MO  +  1TNO+1) )* 
t  NO 4-1) /MO)  *H)-1  .0)  *P1(  2tM0i)  )) 

DGO    TO    200 
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507 


509 


200 
210 
300 


50 


K1G  =B<X,Y+1)*P2(L-1,X) 

K2G=B(X,Yi-l)*0.5*<  P2  (  L  t  X)  *P2(  L-l  ,X  ) ) 

K3G=B(XfY>ll*0.5*(P2(LfX)*-P2(L-ltX)) 

K4G=6(X,Y*1)*P2(L,X) 

P1(L,X)  =  PKL-1  ,  X)MH/6.0)=MklG  +  2.0*K2G*2.0*K3G*K4G) 

GC    TO    200 

K1G  =B(X,Y*1)*P2(L-1,X) 

K2G=8(X,Y*1)*0.5*(P2(L,X)+P2<L-1,X)  ) 

K3G=6<Xf  Y*l)*0.5*(P2<LfX)+P2(L-lfxn 

K4G=bl  A,Yi-l)*P2lL,X) 

P1(L,X)  =  P1(L-1,X)<-(H/6.0)=MK1G  +  2.0*K2G*-2.0*K3G*K4G) 

CCNTINUE 

CCNTINUE 

CONTINUE 

DO    1    I=M3P1,MQ1 

IFdT.EC.OJGO    TO    50 

PRG6U-MBP,Y-NBP)=P1UT,I) 

00    1    J=1,L01 

P2(J,I)=P1( J, I) 

VAR(  J)=VAR(J)*<  I-1.0J*(I-1.0l*PKJfI) 

EST(  J)  =E  SI  (J)  *  (I -1.0)*  PI  (J,  I J 

CONTINUE 

Z=Y*1 

IF(IT.GT.O)G0    TO    901 

IF(M01-MBP.GT.15)G0    TO    901 


C 

c 


c 
c 


901 


99 
400 


2 

850 
851 

700 


110 

900 


WRITE (6, 100) Z, t I ,I=MBP,MO) 


DC    99    J=1,L01 
D( J)=( J-1J*H 


IF(IT„GT.O)GO    TO    99 
IF(MG1-M8P.GT.15 )GG    TO    99 
WRITE(6,1Q1)D( J)  ,(PKJt  I) 


I=M3P1,M01) 


CCNTINUE 

CCNTINUE 

TI  =  (IT-1J.*H 

WPITE(6,105)H 

WPITE(6,103) 

A2=SQRT( A1*B1) 

A3=SQRT(Ai/Bli 

A4=1.0/A3 

A5=M0-NC*A3 

A6  =  MG+i\G*A3 

DC    2    J=1,L01 

LC  =  J 

F1=EXP( A2*D( J) )*A5 

F2=EXP(-A2*D(J) )*A6 

DETERM(J)==0.5*(F1+F2J 

VAR( J)=VA-RU)-cST(  J)*EST(  J) 

WRITE(6,i02)D( J) , EST ( J ) ,DE TERM ( J ) fVAR(JJ 

IFiDETERMI J) .LE.M3P)G0    TO    850 

YF0RCE=0.5*A4*(F2-F1) 

IF(YFGRCE-LE.NBP)GO    TO    851 

CCNTINUE 

GO    TO    700 

Iw  =  l 

GO    TO    700 

IW  =  0 

GO    TO    700 

IF< IT.EQ.OJGO    TO    900 

WRITE( 6f 106JTI 

WRITER  6, 10  7) (NT,NT  =  NBP,NO) 

MJ=MGl-MdP 

NJ=N01-N8P 

DC    110    1=1, MJ 

MT=I+MBP-1 

WRITE(6,104)MT,{ PRCBCI 

IF(IT-GT.OJGU    TO    902 


J) t J=1»NJ) 
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DMIN=DETERM(LC) 

TF1=0.25*D(LC)/H 

TF2=2.0*TF1 

TF3=3.0*TF1 

WRITE (6, 108) M0,N0,MBP,NBP,D<LC) tTFl,TF2,TF3 

1FI  IW.EQ.UGO  TO  852 

WRITE(6,109) 

GO'  TO  8  53 

852  WRITE(6tlll) 

853  IF(EST(LC) .LT. DE TERM( LC) ) GO  TO  270 
DM1=M0-DM4N 

DO  903  J=1,LC 

DETERM(J)  =  U  DETERMi  J J-DMIN ) /DM1 )*100.0 
D< J)=D( J)#100.Q/D(LC) 
903    CGNTINUE 

CALL    PL0T5 

CALL    SCALEiDETERM,  tCl ,9 .0 , 1 .0 , YM  IN, DY ) 

CALL    SCALEID,LC 1,6.0,1.0, XM IN, DX J 

CALL    AXIS(0.0,0.0,XTITLE,-4,6.0,0.0,XHIN,DX) 

CALL    AXIS<0.0,0.0,YTITLE,4,9.0,90.0,YMIN,Dn 

CALL    LINE*D,DETERM,LC, 1,1) 

DO    14    1=1, LC 

ESTU)  =  lUESTl  I)~DMN)/DM1)*100.0-YMIN)/DY 

14  CGNTINUE 

CALL    LINE<D,EST,LC,1,1) 
CALL    PLOTE 
GO    TO    271 

270  DM1=M0-DMIN 
DO    950    J=1,LC 

DETER M( JJ=((DETERM{ J)-DMI N) /DM1) *100 .0 
D( J)=D( J)*10G.0/D(LC) 
950    CGNTINUE 

CALL    PLGT5 

CALL    SCALE(DETERMf  LCI  ,9.0 , 1 .0 , YMI N, DY ) 
CALL    SCALE  ID, LC, 1,6. 0,1.0, XM IN, DX) 
CALL    AXIS{Q.O,0.0,ZTITL£,-4,6.0,G.G,XMIN,DX) 
CALL    AXISIO.O,0.0,YT"ITLE,4,9.0,90.0,0.0,DY) 
CALL    LlNEiD,DETERM,LC,l,i) 
DO    30    1=1, LC 

EST(I)=(( (EST{ I)-DMIN)/DM1)*100.0-YMIN)/OY 
30    CONTINUE 

CALL    LINE(D,EST,LC,1,1) 
CALL    PLOTE 
C 

271  IFUT.EG.OJGO    TO    999 
902    DO    15    K  =  1,.MS 

15  EST(K)=K-1 
DO    16    n=1,NS 

16  DETERM1K)=K-1 
AA=15.0 
BB=30.0 
PMAX=PR0B.(1,  1) 
DC    18    1=1, MS 
DO    18    J=1,NS 

IF(PR03(It J) .GT.PMAX)PMAX=PROB(I ,J) 

18  CONTINUE 
CCNST=8.0VPMAX 
DO    19    1=1, MS 
DO    19    J=1,NS 

PROB(I , J)=CGNST*PR06(I  ,J) 

19  CGNTINUE 

CALL  PLT3D1(EST,MS,DETERM,NS,PR0B,AA,BB,F,TTL,SIZE, 
1WK,IDN,KX,.KY,NKXY,0) 
999  GCNTINUE 

100  FGRMAT( 1X////,3X,' THE  FOLLOWING  ARE  FOR  ',I3,«  SURVI1, 
l'VORS  ON  THE  Y  5  I DE « /// , 2GX , • NO.  OF  X  SURVIVORS1//, 
26X,'TI  ME1-,  1517////) 

101  FGRMAT16X, F7. 3, 15F7. 4) 

102  FORMAT ( 15X.F7.3, 15Xf Fl 1.7, 15Xf F 1 1. 7 , 18X ,F1 1 -7) 

103  FCRMATU5X,'  TIME   • , 1 5X , • E XPECT ED  NO.  OF  SURVIVORS1, 
15X,« DETERMINISTIC  EQUIVALENT1 ,5X ,• VARIANCEt /) 

104  F0RMAT(13X,I7,15F7.4//J 
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105  FOPMATC 1X///,5X, »H  IS  • ,F1 1 .7//// ) 

106  FORMAT! IX, /////♦ 3X , 'MATRIX  GF  PROBABILITIES  OF  NUMBER* 
1«  GF  SURVIVORS  AT  TIME  ',F11.7//) 

107  FQPMATt 19X»15i7//J 

108  FORMAT! lX///t5X» • INITIAL  FORCE  SIZE  FOR  X  IS  «tI3/, 
15X,' INI TIAL  FORCE  SIZE  FOR  Y  IS  •  , 13/ , 5X , ■ BREAKPOI NT • , 
2«  FORCE  LEVEL  FOR  X  IS  • f I  3/ ,5X ,' BREAKPOINT  FORCE  •, 
3'LEVEL  FOR  Y  IS  • ,  1 3/, 5X , ■ BATT LE  WILL  TERMINATE  AT  •, 
4«TIME  • ,F11.7///,5X,« TO  COMPUTE  3-0  GRAPHS  AT  15,     •, 
5'50,  AND  75*  TIMES  INTO  THE  BATTLEtIT  SHOULD  3E«/, 10X, 
6*FCR  25%«  ,.F11.7/jl0X,«FOR  50«"  ,F  1 1 .  7/  ,  10X,  •  FOR  752^1 

109  FCRMAT(5X,.'THE  WINNER  IS  X'////) 
111  FORMAT  (5X-f»  THE  WINNER  IS  Y'////J 

END 

FUNCTION  A(I ,J) 

A  =  0.00d*(J-l  .0) 

RETURN 

END 

FUNCTION  B(I ,JI 

B=0. 004*11-1. 01 

RETURN 

END 
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C  COMPUTER    PROGRAM 

C 

c 

C  SOLUTION    TO    SPRINGALLS    RECURSIVE    SOLUTION    FOR    THE 

C  PROBABILITY    ONE    SIOE    WILL    WIN 

C 

c 

DIMENSION    X( 150,40) ,NS(5) ,RS(8) 

NSU  )=4 

NS<2)=5 

NS(3)=10 

NS(4)=2Q 

NS(5)=40 

RS(1)=0.25 

RS(2)=0.5 

RS{^)=1.0 

RS(4)=2.J 

RS(5)=3.0 

RS(6)=4.0 

RS(7)=5.0 

RS<8)=6.0 

FX=1.2 

FY=1.2 

DO    B    Ki  =  l-,5 

FX=FX-0.2 

FY=FY-0.2 

DO    6    L=l,5 

NC=NS(L) 

DO    7    K=l,8 

R=R  S ( K ) 

SM=2.0*NO.*NO*<  l,0-U.0-FY)*(1.0-FY))*R/(  l.O-(l.O-FX) 
1*(  l.O-FX) ) 

MC-SQRT(SM) 

MBP=(l.Q-FX)*MO 
'  NBP=(  i.0-FY)*NO 

MBP1=MBP+I 

N6P1=NBP*1 

MDIFF=MO-MBP 

NDIFF=NG-N8P 

MCOUNT=MOIFF-l 

NCCUNT=N0IFF-1 

X(MDIFF,NDIFFJ=1.0/(A(MO,NO)+B(MOfNO) J 

DO    1    J=l,NCOUNT 

NREAL=NO-J 

NSTORE=NDIFF-J 

X(MDIFF,NSTORE)=(B(MOtNREAL+l)*X(MDIFF,NSTORE+l) )/ 
l(A(MO, NREAL) «-B<MOf NREAL) J 

1  CONTINUE 

DO    2    I=l,MCOUNT 

KREAL=MO-I 

MSTORE=HDIFF-I 

X(MSTORE,NDIFF)=<  A (MRE AL  +  1 , NO ) *X (MSTORE  + 1, NO  IFF )  )/ 
1 (A(MREAL,NG) +B( M RE AL , NO) ) 

DO    3    J=1,NCGUNT 

NREAL=NO-J 

NSTORE=NDIFF-J 

X(MSTORE,NSTORE)=( A( MR EAL+1, NREAL )*X( MS TORE* l,NSTORE) 
H-B(MREAL,NREALH )*X(MSTORE ,NSTORE+l) )/ 
2<A(MREAL,NREAL) +  B ( NREAL , NREAL  )  ) 

3  CONTINUE 

2  CONTINUE 
PMWIN=0.0 
PNWIN=0.0 

DO    4    I=1,MDIFF 

PMWIN  =  PMWIN*-B(  H-MBP,NBP1)*X(I  ,  1) 

4  CONTINUE 

DO    5    J=1,NDIFF 
PNWIN=PNWIN+A(MBPl,J+NBP)*X(lf J) 

5  CCNTINUE 
SUM=PMWIN<-PNWIN 

S=MO*MO*(l.-(l.-FX)*< l.-FX) )/<2.*(l.-(l. -FY )*(!.- FY) ) 
l*NO*NU) 
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WRITE<6,100)M0,MBP,N0,NBP,FX,FY 
W R  I  T  E  (  6  , 1  0 1 )  PM W  I N  ,  PiM  W I  N  f  SUM  ,  S 

7  CCNTINUE 
6    CONTINUE 

8  CCNTINUE 

100  FCfcMAT(5X, • M0= • ,  14  ,  5X , • MBP=« ,  I4,5X,«NQ=» , 14 , 5X f • N3P- • f 
11<*,5X, *FX=*  tF6.2t5X,,FY=l ,F6.2/J 

101  FCHMAT(5X,  •  PMH  I  N=  •  ,F  1 1 .  7  ,  5X  ,  »  PNW  IN-"  ,F  1 1 .  7  ,  5X,  85UM=«  , 
1F11.7,5X,IR=I,F11.7///J 

END 

FUNCTION  All ,J) 

A=C.C08*J 

RETURN 

ENO 

FUNCTION  B(ItJ) 

B=0. 004*1 

RETURN 

ENO 
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